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The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a 
protoneutron star and a collapse-driven supernova explosion. This process represents a path alter- 
native to thermonuclear disruption of accreting white dwarfs in Type la supernovae. In the AIC 
scenario, the supernova explosion energy is expected to be small and the resulting transient short- 
lived, making it hard to detect by electromagnetic means alone. Neutrino and gravitational-wave 
(GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated 
by the need for systematic predictions of the GW signature of AIC, we present results from an ex- 
tensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation 
of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 
progenitor models in axisymmetric rotational equilibrium, with a wide range of rotational configu- 
rations, temperatures and central densities, and resulting white dwarf masses, we extend previous 
Newtonian studies and find that the GW signal has a generic shape akin to what is known as a 
"Type III" signal in the literature. Despite this reduction to a single type of waveform, we show 
that the emitted GWs carry information that can be used to constrain the progenitor and the post- 
bounce rotation. We discuss the detectability of the emitted GWs, showing that the signal-to-noise 
ratio for current or next-generation interferometer detectors could be high enough to detect such 
events in our Galaxy. Furthermore, we contrast the GW signals of AIC and rotating massive star 
iron core collapse and find that they can be distinguished, but only if the distance to the source is 
known and a detailed reconstruction of the GW time series from detector data is possible. Some of 
our AIC models form massive quasi-Keplerian accretion disks after bounce. The disk mass is very 
sensitive to progenitor mass and angular momentum distribution. In rapidly differentially rotating 
models whose precollapse masses are significantly larger than the Chandrasekhar mass, the resulting 
disk mass can be as large as ~ 0.8Mq. Slowly and/or uniformly rotating models that are limited 
to masses near the Chandrasekhar mass produce much smaller disks or no disk at all. Finally, we 
find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation 
to develop a gravito-rotational bar-mode instability. Moreover, many of our models exhibit suffi- 
ciently rapid and differential rotation to become subject to recently discovered low-Srot/jW^I-type 
dynamical instabilities. 

PACS numbers: 04.25.D-, 04.30.Db, 97.60.Bw, 02.70.Bf, 02.70.Hm 



I. INTRODUCTION 

Single stars with main-sequence masses M < 100 
end their nuclear-burning lives as electron-degenerate ob- 
jects or with central electron-degenerate cores. More 
specifically, the end state is a carbon-oxygen or oxygen- 
neon white dwarf (WD) in the case of low-mass stars 
{i.e., with M <6 — 8 Mq), or a degenerate oxygen-neon 
or iron core embedded in an extended non-degenerate 
stellar envelope in the case of more massive stars {i.e., 6 — 
8Mq < M < 100 Mq) (see, e.g., ^MM and references 
therein). Electron-degenerate spherically-symmetric ob- 
jects become unstable to radial contraction once their 
mass exceeds the Chandrasekhar mass which, assuming 
zero temperature and no rotation, is given by Mch = 



1.4575 (ye/0.5)^ Mq , where Ye is the number of electrons 
per baryon, or "electron fraction" [1, The effective 
Chandrasekhar mass Mch,off of a WD or a stellar core in- 
creases somewhat with WD/core entropy {e.g., [1]) and 
can grow considerably by rotation, in which case it is 
limited only by the onset of nonaxisymmetric instability 
{e.g., Sil). 

The iron core of a massive star is pushed over its Chan- 
drasekhar limit by the ashes of silicon shell burning and 
undergoes collapse to a protoneutron star (PNS), accel- 
erated by photodisintegration of heavy nuclei and elec- 
tron capture 0]. In a high-density sub- Mch oxygen- neon 
core of a less massive star, electron capture may decrease 
Mch,off, also leading to collapse [1, Ho]- In both cases, if 
an explosion results, the observational display is associ- 
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atcd with a Type Il/Ibc supernova (SN). 

On the other hand, a carbon-oxygen WD can be 
pushed over its stabihty hmit through merger with or ac- 
cretion from another WD (double-degenerate scenario) 
or by accretion from a non-degenerate companion star 
(single-degenerate scenario). Here, the WD generally ex- 
periences carbon ignition and thermonuclear runaway, 
leading to a Type la SN and leaving no compact rem- 
nant [III]. However, at least theoretically, it is possible 
that massive oxygen-neon WDs^ formed by accretion or 
merger, and, depending on initial mass, temperature, and 
accretion rate, also carbon-oxygen WDs, may grow to 
reach their Mch.eff or reach central densities sufficiently 
high (> 10^'^ — lO^^'gcm"'^) for rapid electron capture 
to take place, triggering collapse to a PNS rather than 
thermonuclear explosion 0, [IS [H, [H, [13, M, M, M, M, 
[22l . [23l . [2^ . [2^. This may result in a peculiar, in most 
cases probably sub-energetic, low-nickel-yield and short- 
hved transient [13, IH, [23, HI, [H, [s^l • This alternative 
to the Type la SN scenario is called "accretion-induced 
collapse" (AIC) and will be the focus of this paper. 

The details of the progenitor WD structure and forma- 
tion and the fraction of all WDs that evolve to AICs are 
presently uncertain. Binary population synthesis mod- 
els [2^ [m, [13] and constraints on r-process nucleosyn- 
thetic yields from previous AIC simulations [13, [11] pre- 
dict AIC to occur in the Milky Way at a frequency of 
~ 10"^ to ~ 10~*yr"^ which is ~ 20 - 50 times less 
frequent than the expected rate of standard Type la 
SNe {e.g., [13, [H, [H, [13])- In part as a consequence 
of their rarity, but probably also due to their short dura- 
tion and potentially weak electromagnetic display, AIC 
events have not been directly observed (but see [3^ , who 
discovered a peculiar type-lb SN in NGC 1032 that can 
be interpreted as resulting from an AIC). 

The chances of seeing a rare galactic AIC are dra- 
matically boosted by the possibility of guiding elec- 
tromagnetic observations by the detection of neutrinos 
and gravitational waves (GWs) emitted during the AIC 
process and a subsequent SN explosion. GWs, simi- 
lar to neutrinos, are extremely difficult to observe, but 
can carry "live" dynamical information from deep inside 
electromagnetically-opaque regions. The inherent multi- 
D nature of GWs (they are lowest-order quadrupole 
waves) makes them ideal messengers for probing multi- 
D dynamics such as rotation, turbulence, or NS pulsa- 
tions [H, [13, [13]. The detection prospects for a GW 
burst from an AIC are significantly enhanced if theoret- 
ical knowledge of the expected GW signature of such an 
event is provided by computational modelling. In re- 



^ Previously, such WDs were expected to have a significant cen- 
tral ■^^Mg mass fraction, hence were referred to as oxygcn- 
neon-magnesium WDs. Recent work based on up-do-date in- 
put physics and modern stellar evolution codes suggests that the 
mass fraction of ^''Mg is much smaller than previously thought 
(e.g., ^MM)- 



verse, once a detection is made, detailed model predic- 
tions will make it possible to extract physical information 
on the AIC dynamics and the properties of the progen- 
itor WD and, hence, will allow "parameter-estimation" 
of the source. 



Early 



spherically-symmetric (ID) simulations of 



42, 43] and more recent axisymmetric (2D) 
28l . [29[ have demonstrated that the dynam- 
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ics of AIC is quite similar to standard massive star core 
collapse: During collapse, the WD separates into a sub- 
sonically and homologously collapsing (v oc r) inner core 
and a supersonically collapsing outer core. Collapse is 
halted by the stiffening of the equation of state (EOS) at 
densities near nuclear matter density and the inner core 
rebounds into the still infalling outer core. An unshocked 
low-entropy PNS of inner-core material is formed. At its 
edge, a bounce shock is launched and initially propagates 
rapidly outward in mass and radius, but loses energy to 
the dissociation of heavy nuclei as well as to neutrinos 
that stream out from the optically thin postshock re- 
gion. The shock stalls and, in the AIC case (but also 
in the case of the oxygen-neon core collapse in super- 
AGB stars [3 HH), is successfully revived by the de- 
position of energy by neutrinos in the postshock region 
(i.e., the "delayed-neutrino mechanism" [1, \^) or by 
a combination of neutrino energy deposition and rnag- 
netorotational effects in very rapidly rotating WDs [29^ . 
But even without shock revival, explosion would occur 
when the WD surface layer is eventually accreted through 
the shock. Following the onset of explosion, a strong 
long-lasting neutrino-driven wind blows off the PNS sur- 
face, adding to the total explosion energy and estab- 
lishing favorable conditions for r-process nucleosynthe- 
sis [11 [13, [H, [!§, [13]. If the progenitor WD was rotat- 
ing rapidly (and had a rotationally-enhanced Mch,off), a 
quasi-Kcplcrian accretion disk of outer-core material may 
be left after the explosion [H]. Metzger et al. [13, [H] re- 
cently proposed that this may lead to nickel-rich outflows 
that could significantly enhance the AIC observational 
display. 

Rotating iron core collapse and bounce is the most ex- 
tensively studied and best understood GW emission pro- 
cess in the massive star collapse context (see, e.g., [4§] 
and the historical overview in [4l|V However, most mas- 
sive stars (perhaps up to ^ 99% in the local universe) 
are likely to be rather slow rotators that develop little 
asphericity during collapse and in the early postbounce 
phase [10, [13, [13 and produce PNSs that cool and con- 
tract to neutron stars with periods above ~ 10 ms and 
parameter /? = i?rot/]M^] < 0.1% [11], where -Brot is the 
rotational kinetic energy and \W\ is the gravitational 
binding energy. This does not only reduce the over- 
all relevance of this emission process, but also dimin- 
ishes the chances for postbounce gravito-rotational non- 
axisymmetric deformation of the PNS which could boost 
the overall GW emission [4l| . Axisymmetric rapidly ro- 
tating stars become unstable to nonaxisymmetric defor- 
mations if a nonaxisymmetric configuration with a lower 
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total energy exists at a given P (see [53| for a review). 
The classical high-/? instability develops in Newtonian 
stars on a dynamical timescale at /? > — 27% 

(the general-relativistic value is /3 > 25% [Hj, [1^). A 
"secular" instability, driven by fluid viscosity or GW 
backreaction, can develop already at /3 > /Jsec — 14% 
[5^. Slower, but strongly differentially rotating stars 
may also be subject to a nonaxisymmetric dynamical in- 
stability at P as small as ~ 1%. This instability at low 
P was observed in a number of recent 3D simulations 
{e.g., m, M, S M, M, M, M, M, Isl), and may be 
related to corotation instabilities in disks, but its nature 
and the precise conditions for its onset are presently not 
understood [H, [HBl- 

Stellar evolution theory and pulsar birth spin estimates 
suggest that most massive stars are rotating rather slowly 
{e.g., [13, [131, but also [111, [b^] for exceptions). Hence, 
rotating collapse and bounce and nonaxisymmetric rota- 
tional instabilities are unlikely to be the dominant GW 
emission mechanisms in most massive star collapse events 
[4ll |. The situation may be radically different in AIC: In- 
dependent of the details of their formation scenario, AIC 
progenitors are expected to accrete significant amounts 
of mass and angular momentum in their pre-AIC evo- 
lutions 0, [m, m, [13, [ell. They may reach values of P 
of up to ~ 10% prior to collapse, according to the re- 
cent work of Yoon & Langer [3, [l^l, who studied the 
precollapse stellar structure and rotational configuration 
of WDs with sequences of 2D rotational equilibria. De- 
pending on the distribution of angular momentum in the 
WD, rotational effects may significantly affect the col- 
lapse and bounce dynamics and lead to a large time- 
varying quadrupole moment of the inner core, resulting 
in a strong burst of GWs emitted at core bounce. In ad- 
dition, the postbouncc PNS may be subject to the high-/? 
rotational instability (see [69l . [70| for an investigation via 
equilibrium sequences of PNSs formed in AIC) or to the 
recently discovered low-/3 instability. 

Most previous (radiation-)hydrodynamic studies of 
AIC have either been limited to ID [1^, [13, [H or 
were 2D, but did not use consistent 2D progenitor mod- 
els in rotational equilibrium [27l |. Fryer, Hughes and 
Holz [7l| presented the first estimates for the GW sig- 
nal emitted by AIC based on one model of [13|- Draw- 
ing from the Yoon & Langer AIC progenitors [3, [23j . 
Dessart et al. [13, [13] have recently performed 2D 
Newtonian AIC simulations with the multi-group flux- 
limited diffusion (MGFLD) neutrino radiation-(M)HD 
code VULCAN/2D [zl [zS, [zi- They chose two rep- 
resentative WD configurations for slow and rapid rota- 
tion with central densities of 5 x 10^° gcm~^ and total 
masses of 1.46Mq and 1.92^/©. Both models were set up 
with the differential rotation law of [3, HI. The 1.46Mq 
model had zero rotation in the inner core and rapid 
outer core rotation while the 1.92Mq was rapidly ro- 
tating throughout (ratio ilmax,initial/^^conter, initial ~ 1-5). 

Dessart et al. [13, [13| found that rapid electron capture 
in the central regions of both models led to collapse to 



a PNS within only a few tens of milliseconds and re- 
ported successful neutrino-driven [28j and magnetorota- 
tional explosions [13| with final values of P {i.e., a few 
hundred milliseconds after core bounce) of ~ 6% and 
26%, for the 1.46M©and 1.92Mo models, respec- 
tively2. The analysis in^, \M, ^ of the GW signal of 
the Dessart et al. models showed that the morphology of 
the AIC rotating collapse and bounce gravitational wave- 
form is reminiscent of the so-called Type HI signal first 
discussed by Zwerger & Miiller [75[ and associated with 
small inner core masses and a large pressure reduction at 
the onset of collapse in the latter's polytropic models. 

In this paper, we follow a different approach than 
Dessart et al. [13, [13|- We omit their detailed and 
computationally-expensive treatment of neutrino radia- 
tion transport in favor of a simple, yet effective delep- 
tonization scheme for the collapse phase [T^I- This sim- 
plification, while limiting the accuracy of our models at 
postbounce times > 5 — 10 ms, (i) enables us to study 
a very large set of precollapse WD configurations and 
their resulting AIC dynamics and GW signals and, im- 
portantly, (ii) allows us to perform these AIC simulations 
in general relativity, which is a crucial ingredient for the 
accurate modeling of dynamics in regions of strong grav- 
ity inside and near the PNS. Furthermore, as demon- 
strated by [13, [73, [z31 , general relativity is required for 
qualitatively and quantitatively correct predictions of the 
GW signal of rotating core collapse. 

We focus on the collapse and immediate postbounce 
phase of AIC and perform an extensive set of 114 2D 
general-relativistic hydrodynamics simulations. We an- 
alyze systematically the AIC dynamics and the proper- 
ties of the resulting GW signal. We explore the depen- 
dence of nonrotating and rotating AIC on the precol- 
lapse WD rotational setup, central density, core temper- 
ature, and core deleptonization, and study the resulting 
PNS's susceptibility to rotational nonaxisymmetric de- 
formation. Furthermore, motivated by the recent work 
of Mctzger et al. [13, [131 , who discussed the possible en- 
hancement of the AIC observational signature by out- 
flows from PNS accretion disks, we study the dependence 
of disk mass and morphology on WD progenitor charac- 
teristics and rotational setup. 

We employ the gene ral-relativistic hydrodynamics 
code CoCoNuT [73,lzi and neglect MHD effects since 
they were shown to be small in the considered phases un- 
less the precollapse magnetic field strength is extremely 
large {B > 10^^ G, e.g., [13, 113, IMI)- We employ a 
finite-temperature microphysical nuclear EOS in com- 
bination with the aforementioned deleptonization treat- 
ment of [76j . The precollapse 2D rotational-equilibrium 
WDs are generated according to the prescription of Yoon 



^ These numbers are for the non-MHD simulations of [Ssj . In the 
MHD models of [29l . an r2-dynamo builds up toroidal magnetic 
field, reducing the overall rotational energy and /3. 
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& Langer 0,[2|. 

The plan of the paper is as follows. In Sec. |TT1 we 
introduce the numerical methods employed and discuss 
the generation of our 2D rotational-equilibrium precol- 
lapse WD models as well as the parameter space of WD 
structure and rotational configuration investigated. In 
Sec. mil wc discuss the overall AIC dynamics and the 
properties of the quasi-Keplerian accretion disks seen in 
many models. Sec. |IV] is devoted to a detailed analy- 
sis of the GW signal from rotating AIC. There, we also 
assess the detectability by current and future GW obser- 
vatories and carry out a comparison of the GW signals 
of AIC and massive star iron core collapse. In Sec. |Vl 
we study the postbouncc rotational configurations of the 
PNSs in our models and assess the possibility for nonax- 
isymmetric rotational instabilities. In Sec. lVIi we present 
a critical summary and outlook. 



II. METHODS AND INITIAL MODELS 
A. The General-Relativistic Hydrodynamics Code 

We perform our simulations in 2 + 1 dimensions us- 
ing the CoCoNuT code [zl \^ which adopts the 
conformally-flat approximation of general relativity [s^. 
This has been shown to be a very good approximation 
in the context of stellar collapse to PNSs [6l, [H, [sl ]. 
CoCoNuT solves the metric equations as formulated in 
[sst using spectral methods as described in [t^. The 
relativistic hydrodynamics equations are solved via a 
finite-volume approach, piecewise parabolic reconstruc- 
tion, and the HLLE approximate Riemann solver [s^. 
CoCoNuT uses Eulerian spherical coordinates {r, 6} 
and for our purposes assumes axisynimetry. For the com- 
putational grid, we choose 250 logarithmically-spaced, 
centrally-condensed radial zones with a central resolu- 
tion of 250 m and 45 equidistant angular zones covering 
90°. We have performed test calculations with different 
grid resolutions to ascertain that the grid setup speci- 
fied above is appropriate for our simulations. The space 
between the surface of the star and outer boundary of 
the finite difference grid is filled with an artificial atmo- 
sphere. We assume a constant density and stationary 
atmosphere in all zones where density drops below a pre- 
scribed threshold of 7 x lO^gcm"'^, a value marginally 
larger than the lowest density value in the EOS table 
employed in our calculations (cf. Sec. Ill A The at- 
mosphere is reset after each timestep in order to ensure 
that it adapts to the time-dependent shape of the stel- 
lar surface. For further details of the formulations of 
the hydrodynamics and metric equations as well as their 
numerical implementation in CoCoNuT, the reader is 
referred to [4i,[7i,[8|. 

The version of CoCoNuT employed in this study does 
not include a nuclear reaction network. Hence, we, like 
Dessart ct al. (28l.[29j. ignore nuclear burning which may 
be relevant in the outer core of AIC progenitors where 



material is not in nuclear statistical equilibrium (NSE), 
but still sufficiently hot for oxygen/neon/magnesium 
burning to occur. This approximation is justified by re- 
sults from previous work of (26l . |4^ who included nuclear 
burning and did not observe a strong dynamical effect. 



1. Equation of State 

We make use of the finite-temperature nuclear EOS of 
Shen et al. ("Shen-et-al EOS" in the following, [13, IH) 
which is based on a relativistic mean-field model and is 
extended with the Thomas-Fermi approximation to de- 
scribe the homogeneous phase of matter as well as the 
inhomogeneous matter composition. The parameter for 
the incompressibility of nuclear matter is 281 MeV and 
the symmetry energy has a value of 36.9 MeV. The Shen- 
et-al EOS is used in tabulated fashion and in our version 
(equivalent to that used in [i^, [s^) includes contribu- 
tions from baryons, electrons, positrons and photons. 

The Shen-et-al EOS table used in our simulation 
has 180, 120, and 50 equidistant points in logj^Q p, 
log^Q T, and Ye, respectively. The table ranges are 
6.4x lO'^gcm-^ <p< l.lx lO^Sgcm-^, O.lMcV < T < 
100.0 MeV, and 0.015 < < 0.56. Our variant of the 
Shen-et-al EOS assumes that NSE holds throughout the 
entire {p,T,Yf,} domain. In reality, NSE generally holds 
only at r > 0.5 MeV. At lower temperatures, a nuclear 
reaction network and, the advection of multiple chemi- 
cal species and accounting for their individual ideal-gas 
contributions to the EOS is necessary for a correct ther- 
modynamic description of the baryonic component of the 
fluid. However, since the electron component of the EOS 
is vastly dominant in the central regions of AIC progeni- 
tors (and also in the central regions of iron cores), the in- 
correct assumption of NSE at low temperatures can lead 
to only a small error in the overall (thermo)dynamics of 
the collapse and early postbounce phase. 



2. Deleptonization during Collapse and Neutrino Pressure 

To account for the dynamically highly important 
change of the electron fraction by electron capture 
during collapse, we employ the approximate prescription 
proposed by Liebendorfer (76| . Liebendorfer's scheme is 
based on the observation that the local Ye of each fluid 
element during the contraction phase can be rather ac- 
curately parametrized from full radiation-hydrodynamics 
simulations as a function of density alone. Liebendorfer 
demonstrated the effectiveness of this parametrization 
in the case of spherical symmetry, but also argued that 
it should still be reliable to employ a parametrization 
Ye{p) obtained from a ID radiation-hydrodynamics cal- 
culation in a 2D or 3D simulation, since electron cap- 
ture depends more on local matter properties and less 
on the global dynamics of the collapsing core. On the 
basis of this argument, a Ye{p) parametrization was ap- 
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FIG. 1: Average electron fraction in the innermost 2 km in 
the collapsing WD as a function of density obtained from 2D 
MGFLD simulations with the VULCAN/2D code for models 
1.46Mo and 1.92M0 of Dessart et al. Q. Both models were 
set up with the same initial dependence of temperature on 
density and a temperature To = 1.0 x 10^" K (see Sec. Ill Bl 
for details). 



jcctorics obtained from VULCAN/2D AIC simulations. 
At nuclear density, these data predict Yg ^ 0.18 which is 
low compared to y. > 0.22 — 0.26 seen in simulations of 
iron core collapse [U [tI, [qO, [11] and oxygen-neon core 
collapse Q. This difference is not fully understood, but 
(j) could be physical and due to the WD initial data used 
here and in [28| or (ii) may be related to the radiation 
transport approximations and microphysics treatment in 
VULCAN/2D. To measure the importance of these un- 
certainties in YJIp), we perform calculations with system- 
atic variations of Ye{p) due cither to changes in the pre- 
collapse WD temperature or to an ad-hoc scaling (see 
Sec. lirC3l) . 

Since AIC progenitors may be extremely rapidly ro- 
tating, it is not clear that the Ye(p) parametrization is 
indeed independent of the specific model and rotational 
setup. The Ye{p) trajectories shown in Fig. [T] result from 
the collapse simulations of the slowly rotating 1.46 A/0 
model and of the rapidly rotating 1.92 Mq model of 
Dessart et al. The very close agreement of the two curves 
suggests that rotational effects have only a small influ- 
ence on the prebounce deleptonization and confirm the 
supposition of [tHI at the level of the MGFLD and micro- 
physics treatment in VULCAN/2D. All Ye{p) data used 
in this study are available from [o^ ]. 



plied in the rotating iron core collapse calculations of 
pgl . [6^ . m, [il, [z!]. Here, we use the same implemen- 
tation as discussed in [4§| and track the changes in Ye 
up to the point of core bounce. After bounce, we simply 
advect Y^. Furthermore, as in [ioj . we approximate the 
pressure contribution due to neutrinos in the optically- 
thick regime (p > 2 x 10^^ g cm~^) by an ideal Fermi gas, 
following the prescription of [t^. This pressure contri- 
bution and the energy of the trapped neutrino radiation 
field are included in the matter stress-energy tensor and 
coupled with the hydrodynamics equations via the energy 
and momentum source terms specified in [63l |. 

The deleptonization scheme described here is applica- 
ble only until core bounce and can neither track the post- 
bounce neutrino burst (see, e.g., [l^l) nor neutrino cool- 
ing/heating and the postbounce deleptonization of the 
PNS. The dynamics in the very early postbounce evo- 
lution (up to ^ 10 ms) are unlikely to be dramatically 
affected by this limitation, but it should be kept in mind 
when interpreting results from later postbounce times. 

For our AIC simulations we obtain Ye{p) data from 
the 2D Newtonian radiation-hydrodynamics simulations 
carried out by Dessart et al. [2^ with the VULCAN/2D 
code [zi. III, [in in its MGFLD variant. We use these 
data because of their ready availability, but point out 
that the microphysics [9l| used in VULCAN/2D does 
not yet include the updated electron capture rates of 
[9^ . Moreover, VULCAN/2D presently does not treat 
velocity-dependent terms in the transport equation and 
neglects neutrino-electron scattering, both of which my 
have some impact on the evolution of Ye in the collapse 
phase [9O, [9I]. In Fig. [T] we plot representative Ye{p) tra- 



B. Precollapse White Dwarf Models 

For constructing 2D WD models in rotational equilib- 
rium with a given rotation law, we follow ^23*1 and employ 
the self-consistent field (SCF) method ^, M., M, M in 
Newtonian gravity. For the purpose of the SFC method, 
we assume that the WD is cold and has a constant Ye 
of 0.5. After finding the 2D equilibrium configuration, 
we impose a temperature and Ye distribution motivated 
by previous work [H, [2^ . Ideally, the WD initial model 
should be evolved in a multi-dimensional stellar evolu- 
tion code with a finite-temperature EoS and accounting 
for weak processes such as neutrino cooling and electron 
capture. Due to the unavailability of such self-consistent 
AIC progenitors, we resort to the treatment that we dis- 
cuss in detail in the remainder of this section. 



1. Implementation of the Self-Consistent Field Method 

Our implementation of the Newtonian SFC method 
has been tested by reproducing the WD models presented 
in [H, [9I], and finding excellent agreement. The com- 
pactness parameter GM/Rc^ of the highest-density WD 
models considered here reaches 5x 10"'^, hence gencral- 
relativistic effects at the precollapse stage are small and 
the error introduced by Newtonian WD models is there- 
fore negligible. Hereafter we will assume that the Newto- 
nian mass of the equilibrium model represents the haryon 
mass accounted for when solving the general-relativistic 
equations. 
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The equation governing the steUar equiUbrium is given 



by 



(1) 



where $ is the gravitational potential, Q is the angular 
velocity, zu is the radial cylindrical coordinate and C is 
a constant that will be determined from boundary con- 
ditions using the SCF iterations as discussed below. 

White dwarfs are stabilized against gravity by electron 
degeneracy pressure. For constructing precollapse WDs, 
we assume complete degeneracy for which the WD EOS 
{e.g., dl) is given by 



P = A[a;(2x2-3)(x2 + l)i/^ + 3sinh-^x] 



(2) 

6.01 X 10^2 dyn cm^^ ^^^^ ^ ^ g g2 x 
We set Ye = 0.5, assuming at this stage 
that no electron capture has taken place. The integral 
J p^^dP in Eq. ([T]) is the enthalpy H which, given our 
choice of WD EOS, can be expressed analytically as 



where A 

10^ Y-^ g cm-3 



H = 



8A 



2/3' 



With this, Eq. (jT]) trivially becomes 



7? ==C-$ 



n VJ dvj 



(3) 



(4) 



Following the SCF method, we proceed to first produce 
a trial density distribution p(r, 6) and impose a rotation 
law (discussed in the following Sec. Ill B"2|) . 

We then calculate C by using the value for the maxi- 
mum density and the angular velocity at the center of the 
star n{w = 0) = f^ci- Based on the trial density distri- 
bution, we calculate H via Eq. Q and then update the 
density distribution based on H using the analytic ex- 
pression ([2]). This updated density distribution in turn 
results in a new value for H. We iterate this procedure 
until all the maximum absolute values of three relative 
differences of H, fl and p become less than 10"'^. 



2. Progenitor Rotational Configuration 

Our axisymmetric progenitor WD models are assumed 
to be either in uniform rotation or to follow the differen- 
tial rotation law proposed by Yoon & Langer |23| . The 
latter argued that the rotation law of a WD that ac- 
cretes matter at high rates (> 10~^ Mq yr"-'^) is strongly 
affected by angular momentum transport via the dynam- 
ical shear instability (DSI) in the inner region, and due 
to the secular shear instability (as well as Eddington- 
Sweet circulations in the outer layers. According 
to their results, the shear rate in the core remains near 
the threshold value for the onset of the DSI. This results 
in a characteristic rotation law which has an absolute 



maximum in the angular velocity just above the shear- 
unstable core. We define zup as the position of this max- 
imum. This position is linked to layers with a density as 
low as several percent of the WD central density so that 



Piizu 



G7o 



Z ^0) = fpPc 



(5) 



and where, following [23,l2^, we choose /p = {0.05, 0.1} 
in our models. (Note that the differential rotation law 
adopted for the models of [111 had /p = 0.05). In the 
inner regions with zu < Wp, we have 



/shCDSI^crit 



dvo', 



(6) 



where 51c, i is the angular velocity at the center and 
CDSi,ciit is the threshold value of the shear rate in the 
inner core for the onset of the DSI. fsh is a dimensionless 
parameter (< 1.0) describing the deviation of the shear 
rate from (iDSi.crit- We compute CDSi.crit assuming homo- 
geneous chemical composition and constant temperature, 
in which case ctdsi crit can be estimated as (c/. Eq.(7) 
of 0): 



0.2 
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109 cm s- 



(7) 



0.01 



X 10"^ cm 



Vad 

0.4 



where g is the free-fall acceleration. Hp is the pressure 
scale height (= —dr/d\nP), Vad is the adiabatic tem- 
perature gradient (= — {dlnT/ din P)s where s is the 
specific entropy) and S = (dhip/dhiT)p. The quan- 
tities 6, Hp and Vad are computed using the routines of 
Bhnnikov et al. 

At the equatorial surface, the WD is assumed to rotate 
at a certain fraction /k of the local Keplerian angular 
velocity Ok: 



(8) 



where Re is the equatorial radius of the WD and where 
we have set /k = 0.95. In the region between Wp and 
i?e, we again follow [l^ and adopt the following rotation 
law: 

n{m)/nK = n{-ujp) / nKi^Up) +C{vu - rup)" , (9) 

where the constant C is determined for a given value of 
a as 



C = 



/k — r^(tx7p)/r^K(g7p) 

{Re - Wp)° 



(10) 



The choice of the exponent a does not have a strong 
impact on the WD structure because of the constraints 
imposed by fl{-cup) and fl{Re) at each boundary. In our 
study, we adopt a = 1.2. For further details, we refer the 
reader to Sec. 2.2 o f |2 3l. 

Saio & Nomoto [22| argued that turbulent viscosity 
resulting from a combination of a baroclinic instability 
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(see, e.g., jlO0j |: neglected by Yoon & Langer 0, H^) 
and the DSI is so efficient in transporting angular mo- 
mentum that the angular velocity becomes nearly uni- 
form in the WD interior, while only surface l ayers with 
mass < O.OIMq rotate differentially Piro [loH . who 
also considered angular momentum transport by mag- 
netic stresses, confirmed these results. Hence, in order 
to study the suggested case of uniform precollapse WD 
rotation, we complement our differentially rotating WD 
models with a set of uniformly rotating AIC progenitors. 

3. Initial Temperature Profile 

Because our initial models are constructed by impos- 
ing hydrostatic equilibrium (Eq. ([T])) with a barotropic 
EOS (Eq. ([2])), the WD structure is independent of tem- 
perature. However, the latter is needed as input for the 
finite-temperature nuclear EOS used in our AIC simula- 
tions. We follow Dessart et al. [1^ and impose a scaling 
of the temperature with density according to 

T(n7, z) = To [pc,i/p(^, ^)]°'^% (11) 

where (txj, z) are cylindrical coordinates and po is the 
density at which the stellar temperature equals Tq. 

4. Initial Electron Fraction Profile 

For the purpose of constructing AIC progenitor WDs 
in rotational equilibrium we assume that no electron cap- 
ture has yet taken place and set Ye = 0.5. A real AIC 
progenitor, however, will have seen some electron cap- 
tures on Ne/Mg/Na nuclei {e.g., [l^l) before the onset 
of dynamical collapse. In addition, electrons will be cap- 
tured easily by free protons that are abundant at the 
temperatures of the WD models considered here. Hence, 
a Ye of 0.5 is rather inconsistent with real WD evolu- 
tion. Dessart et al. [2^, who started their simulations 
with Ye = 0.5 models, observed an early burst of elec- 
tron capture. This led to a significant initial drop of 
Ye that leveled off after 5 — 10 ms beyond which the Ye 
profile evolved in qualitatively similar fashion to what is 
known from iron core collapse (see Fig. [1] which depicts 
this drop of Ye at low densities). To account for this, we 
adopt as initial Ye{p) a parametrization obtained from 
the equatorial plane of the models of Dessart et al. [Ill 
at ~ 7 ms into their evolution when the initial electron 
capture burst has subsided. We use these Ye{p) data for 
the Ye evolution of the low-density (p < pc.i) part of the 
WD during collapse. 

C. Parameter Space 

The structure and thermodynamics of the AIC progen- 
itor and the resulting dynamics of the collapse depend on 
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FIG. 2: The radial profile of the rest-mass density for nonro- 
tating white dwarf models AUO, BUO, CUO and DUO. 

a variety of parameters that are constrained only weakly 
by theory and observation {e.g., 0, [H, HI]). Here, we 
study the dependence on the central density, rotational 
configuration and the core temperature. In the follow- 
ing we lay out our parameter choices and discuss the 
nomenclature of our initial models whose key properties 
we summarize in Table HI 

1. Progenitor White Dwarf Central Density 

In order to investigate the impact of the precollapse 
central WD density pc,i on the collapse dynamics, we 
consider sequences of WD models with central densities 
in the range from 4 x 10^ gcm^^^ to 5 x 10^" gcm^'^. This 
range of densities is motivated by previous studies argu- 
in g th at WDs in this range of pc i may experience AIC 

We therefore choose a set of four central densities. 
I.e., 4 X 10^ 1 X 10^°, 2 X 10^°, 5 x 10^° gcm'^, 
and correspondingly begin our model names with let- 
ters A, B, C, D. We perform AIC simulations of nonro- 
tating (spherically-symmetric) WDs with central density 
choices A-D and restrict the rotating models to the lim- 
iting central density choices A and D. 

In Fig. [5] we plot radial density profiles of our nonro- 
tating WD models to show the strong dependence of the 
WD compactness on the choice of central density. This 
aspect will prove important for the understanding of the 
collapse dynamics of rapidly rotating models. 

2. Progenitor White Dwarf Rotational Configuration 

Since the rotational configuration of AIC progenitor 
WDs is constrained only poorly, we consider uniformly 
rotating (f2i = rJc.i everywhere) as well as a variety of 
differentially rotating WD configurations. To denote the 
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FIG. 3: Upper panels: angular velocity as a function of equa- 
torial radius (left panel) and enclosed mass coordinate (right 
panel) for three representative precollapse WD models ADS, 
AD5 and ADIO. Lower panels; Angular velocity normalized 
to the local Keplerian value as a function of equatorial radius 
and enclosed mass for the same models. 

general rotation type, we use the letter U (D) for uniform 
(differential) rotation as the second letter in each model 
name. 

The low-density uniformly rotating model sequence 
AU{l-5} is set up with initial angular velocities r^c.i from 
1 to 3.5 rads~^, where the latter value corresponds to ro- 
tation very close to the mass-shedding limit. The more 
compact uniformly-rotating sequence DU{l-7} is set up 
with precollapse Oc,i from 2 to 9.5rads^^, where, again, 
the latter value corresponds to near-mass-shedding rota- 
tion. 

Model sequences AD{1-10}, DD{l-7} are differen- 
tially rotating according to the rotation law discussed in 
Sec. lIIB^ and specified by Eqs. ([6]) and ([9]), with the pa- 
rameter choice /sh — 1 and /p = 0.1 for the AD sequence 
and /sh = 1 and /p — 0.05 for the DD sequence. We re- 
call that /p is the fraction of the central density where the 
angular velocity has a global maximum. While /p = 0.1 
is the standard choice of [1^, we adopt /p = 0.05 for the 
high-density sequence DD to be in line with the param- 
eter choices made for the models of Dessart et al. psl ]. 
Test calculations with AD models show that the varia- 
tion of /p between 0.05 and 0.10 affects the rotational 
configuration of the outer WD layers only and does not 
have any appreciable effect on the AIC dynamics. For 
the AD{1-10} sequence, we chose ric.i in the range from 
to 5.6 rads~^, resulting in maximum angular velocities 
f^max,i in the range of 2.88 to 8.49 rads~^. The higher- 
density DD{l-7} sequence rotates with J7c,i in the range 
from to 18rads~^, corresponding to maximum f2 in 
the range of 7.69 to 25.84 rads~^. The values of Slc,i and 
^^max.i for thc individual AD and DD models are given in 
Table HI As representative examples resulting from our 
assumed rotation law, we plot in Fig. [3] for models ADS, 



AD5, and ADIO the angular velocity and thc ratio of the 
angular velocity to thc local Keplerian value as a func- 
tion of cylindrical radius and of thc enclosed rest mass. 
In Fig. m we plot the colormaps of thc rest mass density 
on the r — Q plane for the representative precollapse WD 
models ADl, AD5 and ADIO. 

In order to study the effect of variations in thc degree 
of differential rotation, we vary the dimensionless shear 
parameter /sh for a subsequence of AD models and ap- 
pend suffices /{I — 4} to their names corresponding to 
/sh = {0.8,0.6,0.4,0.2}, respectively. Figure [5] shows the 
behavior of the initial angular velocity distribution with 
decreasing /sh in the rapidly differentially rotating model 
ADIO. 

An important point to mention is the large range of 
precollapse WD masses covered by our models. Depend- 
ing on thc initial central density and the rotational setup, 
our WDs masses range from a sub-Afch value of 1.39 Mq 
in thc nonrotating low-pc,i model AUO to a rotationally- 
supported super- Mch mass of 2.05 Mq in the rapidly dif- 
ferentially rotating model AD13f4. The maximum mass 
in our sequence of uniformly rotating WDs is 1.462 A/q 
and is obtained in model DU7. 

To conclude the discussion of our initial rotational con- 
figurations, we present in Fig. [6] for all models thc initial 
values (/3i) of thc parameter /3 as a function of their pre- 
collapse central angular velocity fic.i- Differentially ro- 
tating WD models can reach /3i of up to ^ 10% while 
staying below the mass-shedding limit. This number is 
more than a factor of 2 larger than what seems possible 
in massive star iron core collapse (see, e.g., [49l|). making 
these rapidly rotating AIC progenitor models potential 
candidates for a dynamical nonaxisymmetric rotational 
instability during their postbounce AIC evolution (see 
Sec. El). 



3. Progenitor White Dwarf Core Temperature and Ye{p) 
parametrization 

We use Eq. (fTTj) to set up the initial temperature 
distribution as a function of density. Dessart et al. 
chose pq = pc,i (= 5 X 10^" g cm"'^ in their models) and 
To = 10^° K for their 1.46M0 model, and Tq = 1.3 x 10^° 
K for their 1.92AfQ model. These values (i) are simi- 
lar to what was used in the earlier work of Woosley & 
Baron |2^ and (ii) work well with the tabulated EOS em- 
ployed and the assumption of NSE, but may be higher 
than the temperatures prevailing in accreting precollapse 
WDs in nature (see, e.g., 0,[l3,|2l). 

While the fluid pressure is affected very little by dif- 
ferent temperature distributions, this is not the case for 
the free proton fraction which increases strongly with T 
in the range from 10^ to 10^° K and at precollapse core 
densities. This increase of the proton fraction can lead 
to enhanced electron capture during AIC and in this way 
may have a significant influence on thc AIC dynamics. In 
order to test the sensitivity of our AIC simulations on the 
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FIG. 4: Colormap of the rest mass density for the precoUapse white dwarf models ADl (right panel) AD5 (center panel), and 
AD 10 (right panel). The apparent ruggedness of the WD surface layers is a results of the finite resolution of our computational 
grid and the mapping procedure in the visualization tool. The ruggedness has no influence on the collapse and postbounce 
dynamics of the inner core. 




FIG. 5: Angular velocity as a function of equatorial radius for 
model ADIO and varying values of the dimensionless shear 
parameter f^h that controls the rate at which the angular 
velocity increases with vj in the region zu < vjp. 



assumed Tqj we not only study models with Tq = 10^'^ 
K (at po = 5 X lO^^gcm"^, herafter the "high-T" mod- 
els), but perform also simulations for models set up with 
To = 5 X 10^ K (at po = 5 x 10i°_gcm~^ hereafter the 
"low-T" models). To obtain the Yj^p) parametrization 
(sec Sec. Ill for the latter temperature, we re-ran with 
VULCAN/2D the IAQMq AIC model of Dessart at al. 
up to eore bounce with the same setup as discussed 
in [2I], but using the lower value of Tq- We do not indi- 
cate the two different initial temperatures in the model 
names, but list the results obtained in the two cases side- 
by-side in Table mi 



FIG. 6: Parameter /3i versus central angular velocity f2c,i for 
our AIC progenitor WD model set. 



In addition to variations in deleptonization due to dif- 
ferences in the precollapse WD thermodynamics, we must 
also consider the possibility of unknown systematic biases 
that lead to small values of Kg in the inner core at bounce 
(see Sec. Ill A"2|) . In order to study the effect that larger 
values of Yp in the inner core have on the AIC dynam- 
ics, we perform a set of test calculations with scaled Ye{p) 
trajectories. We implement this by making use of the fact 
that Ye{p) is to good approximation a linear function of 
log(p) (see Fig. [1]). We change the slope of this function 
between p = 5 x 10^*'gcm~'^ and p ^ 2.5 x 10^'^gcm^'^ 
by increasing Yeip = 2.5 x lO^^'gcm-^) by 10% and 20%. 
We pick these particular scalings, since the 20% increase 
yields inner-core values of Ye at bounce that are very 
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close to those obtained in ID Boltzmann neutrino trans- 
port simulations of oxygen-neon core collapse [3, Il02ll . 
The 10% scaling yields values in between those of 
and [3, Il02| and, hence, allows us to study trends in 
AIC dynamics with variations in deleptonization in be- 
tween constraints provided by simulations. We will not 
list the results of these tests in our summary tables, but 
discuss them wherever the context requires their consid- 
eration {i.e., Sect. nira hitbi nvi irvTi . 



where Sh is the power spectral density of the detec- 
tor and Sh,c = Shifc)- We approximate the average 
(/i^) over randomly distributed angles by (3/2) h^. From 
Eqs. (fT4|) . (|15|1 the optimal single-detector signal-to-noisc 
ratio (SNR) can be calculated as 

SNR^— (16) 

where /irms = V f^h is the value of the root-mean-square 
strain noise for the detector. 



D. Gravitational Wave Extraction 

We employ the Newtonian quadrupole formula in the 
first-moment of mome ntum density formulation as dis- 
cussed, e.g., in [tqI. [s^ . [Tosj . In essence, we compute the 
quadrupole wave amplitude A^q of the / = 2, m = mode 
in a multipole expansio n of the radiation field into pure- 
spin tensor harmonics |104| |. In axisymmetric AIC, this 
quadrupole mode provides by far the largest contribution 
to the GW emission and other modes arc at least one 
or more orders of magnitude smaller. Of course, should 
nonaxisymmetric instabilities develop (which we cannot 
track in our current 2D models), these would then pro- 
vide a considerable nonaxisymmetric contribution to the 
GW signal. 

The GW amplitude is related to the dimcnsionless GW 
strain h in the equatorial plane by 



1 /15 



AE2 
^20 



= 8.8524x10" 



21 ( \ kpc 

V 103 cm/ V r 

(12) 

where r is the distance to the emitting source. 

We point out that although the quadrupole formula is 
not gauge invariant and is only valid in the weak-field 
slow- motion limit, it yields results that agree very well in 
phase and to ~ 1 0-20% in am plitude with more sophis- 
ticated methods [lol ITosl . [l06t . 

In order to assess the prospects for detection by cur- 
rent and planned interferometric detectors, we calculate 
characteristic quantities for the GW signal following [3^ . 
Performing a Fourier transform of the dimcnsionless GW 
strain h, 



f^hdt. 



(13) 



we can compute the (detector dependent) integrated 
characteristic frequency 




Sh 



-df 



(14) 



and the dimensionless integrated characteristic strain 

1/2 



hr.= 3 



Sh,c /J- 2 







{h')fdf 



(15) 



III. COLLAPSE DYNAMICS 

The AIC starts when the progenitor WD reaches its 
effective Chandrasekhar mass and pressure support is re- 
duced due to electron capture in the core. Similar t o the 
case of massive star iron core collapse {e.g., [tsI. [s^ . [iOTI 
Il08l | and references therein), the collapse evolution can 
be divided into three phases: 

Infall: This is the longest phase of collapse and, de- 
pending on model parameters, lasts between ~ 25 ms and 
^ 300 ms. The inner part of the WD core (the "inner 
core"), which is in sonic contact, contracts homologously 
{vr oc r), while the "outer core" collapses supersonically. 
Fig. [7] shows the time evolution of the central density for 
the nonrotating high-T AIC models. In the infall phase, 
the core contracts slowly, which is reflected in the slow 
increase of pc. 

Plunge and bounce: The short dynamical "plunge" 
phase sets in when pc reaches ~ lO^^gcm"^^ and the 
peak radial infall velocity is ~ 0.1c. At this point, neu- 
trinos begin to be trapped in the inner core. The lat- 
ter rapidly contracts to reach nuclear densities (pnuc ~ 
2.7 X 10^'' g cm~3) at which the nuclear EOS stiffens, 
decelerating and eventually reversing the infall of the 
inner core on a millisecond timescale. Because of its 
large inertia and kinetic energy, the inner core does 
not come to rest immediately. It overshoots its equi- 
librium configuration, then bounces back, launching a 
shock wave at its outer edge into the still infalling outer 
core. The bounce and the re-expansion of the inner 
core is also evident in the time evolution of pc shown 
in Fig. [7| which, at core bounce, reaches a value of 
~ 2.8 X lO^'' gcm~3 in the nonrotating AIC models, after 
which the core slightly re-expands and settles down at 
~ 2.5 X 10^'*gcm^3^ As pointed out by extensive previ- 
ous work (see, e.g., [1, \m ^M, EM 

lllll and references 
therein), the extent of the inner core at bounce deter- 
mines the initial kinetic energy imparted to the bounce 
shock, the mass cut for the material that remains to be 
dissociated, and the amount of angular momentum that 
may become dynamically relevant at bounce. 

Ringdown: Following bounce, the inner core oscillates 
with a superposition of various damped oscillation modes 
with frequencies in the range of 500 — 800 Hz, exhibiting 
weak low-amplitude variations in pc (Fig. [7]). These oscil- 
lations experience rapid damping on a timescale of 10 ms 
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TABLE I: Summary of the initial WD models: Qc,i is the central angular velocity and fimax,! = f2(ti7p), Mq is the total rest-mass 
and J is the total angular momentum. \Wi\ and i7rot,i are the gravitational energy and rotational kinetic energy of the WD, 
respectively. Re and -Rp are the equatorial and polar radii. 
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0.47 
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AD3f3 
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0.26 


38.56 


0.36 


0.90 


2298 
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AD3f4 
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1.883 


0.4 


1.417 


0.21 


38.20 


0.26 


0.70 


2243 


0.745 
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0.4 
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0.64 
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1.472 
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1.89 


2462 


0.655 


AD9fl 


5.227 


7.564 


0.4 


1.772 


1.23 


50.92 


3.71 


7.28 


3574 


0.400 


AD9f2 


5.227 


6.978 


0.4 


1.691 


1.00 


48.13 


2.89 


6.01 


3264 


0.448 


AD9f3 


5.227 


6.392 


0.4 


1.630 


0.84 


46.00 


2.29 


4.98 


3029 


0.493 


AD9f4 


5.227 


5.808 


0.4 


1.584 


0.71 


44.38 


1.83 


4.12 


2851 


0.533 


ADlOfl 


5.554 


7.896 


0.4 


1.833 


1.41 


53.06 


4.35 


8.20 


3793 


0.370 


AD10f2 


5.554 


7.305 


0.4 


1.741 


1.13 


49.93 


3.37 


6.74 


3434 


0.420 


AD10f3 


5.554 


7.721 


0.4 


1.665 


0.93 


47.28 


2.64 


5.58 


3149 


0.468 


AD10f4 


5.554 


6.134 


0.4 


1.611 


0.78 


44.40 


2.10 


4.62 


2942 


0.510 


ADllf2 


6.000 


7.756 


0.4 


1.815 


1.35 


52.60 


4.16 


7.90 


3696 


0.380 


AD12f3 


7.000 


8.175 


0.4 


1.914 


1.65 


56.27 


4.23 


9.30 


4010 


0.340 


AD12f4 


7.000 


7.586 


0.4 


1.798 


1.29 


52.29 


3.99 


7.64 


3574 


0.393 
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TABLE II: Summary of key quantitative results from our AIC simulations. pmax,b is the maximum density in the core at the 
time of bounce th, |/i|max is the peak value of the GW signal amplitude, while /3ic,b is the inner core parameter /3 at bounce. 
Models marked by unfilled/filled circles (•/o) undergo a pressure-dominated bounce with/without significant early postbounce 
convection. Models marked with the cross sign (x) undergo centrifugal bounce at subnuclear densities. The values left/right 
of the vertical separator (|) are for the models with low/high temperature profiles. 
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# 
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• 


• 
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due to the emission of strong sound waves into the post- 
shock region which steepen into shocks. The newly born 
PNS thus rings down to its new equihbrium state. 



The ringdown phase is coincident with the burst 
of neutrinos that is emitted when the bounce shock 
breaks out of the energy-dependent neutrinospheres (see, 
e.g., [28l. [90|). The neutrino burst removes energy from 
the postshock regions and enhances the damping of the 
PNS ringdown oscillations {e.g., [H^l), but, due to the 
limitations of our present scheme (see Sec. Ill A"2|) . is not 
accounted for in our models. 



A. Nonrotating AIC 



The set of nonrotating AIC models that we consider 
here consists of mode ls AUO, BUO, CUO, DUO. As 
noted in Section UlC 1[ these models have different cen- 
tral densities with values in the range from 4 x 10^ to 
5 X 10^° g cm^'^ which, because of the strong dependence 
of the WD compactness on the central density, corre- 
sponds to a range of WD radii from 1692 to 813 km 
(see Fig. [2]). Once mapped onto our computational grid 
and after the initial Ye (p) parametrization is applied (see 
Sec. IIIB4[) . all WD models start to collapse by them- 
selves and no additional artificial pressure reduction is 
necessary. This is in contrast to previous work that em- 
ployed a simple analytic EOS and required an explicit 
and global change of the adiabatic exponent to initiate 
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FIG. 7: Evolution of the maximum (central) density for the 
nonrotating low-T models AUO, BUO, CUO and DUO. The 
inset plot displays a zoomed-in view of the ma^cimum den- 
sity around the time of core bounce on a linear scale. As 
clearly discernible from this figure, the collapse dynamics in 
the plunge and bounce phase are essentially independent of 
the initial WD central density. Time is normalized to the time 
of bounce th- 



collapse {e.g., fldi\7l\). 

The free-fall collapse time Tff of a Newtonian self- 
gravitating object of mean density Pmoan is proportional 

— 1/2 

to pmean. For our sct of splierlcally-symmetric AIC mod- 
els we find a scaling Tff oc p~^'^'^ , where pc is the precol- 
lapse central density of the WD. This stronger scaling is 
due to the fact that WD cores are not constant density 
objects and that the collapse is not pressureless. Further- 
more, the pressure reduction initiating and accelerating 
collapse is due primarily to electron capture which scales 
roughly with p^^^ {e.g., §]). Hence, lower-density WDs 
collapse only slowly, spending much of their collapse time 
near their initial equilibrium states. 

In Fig. [71 we plot the evolution of the central densities 
of the nonrotating high-T models. Despite the strong 
dependence of the collapse times on the initial central 
densities, the evolution of pc around bounce does not ex- 
hibit a dependence on the initial central density. More- 
over, the mass and the size of the inner core is rather 
insensitive to the initial value of pc. 

These features, somewhat surprising in the light of the 
strong dependence of the collapse times on the initial 
value of Pc, are a consequence of the fact that the inner 
core mass is determined primarily not by hydrodynamics, 
but by the thermodynamic and compositional structure 
of the inner core set by nuclear and neutrino physics Q. 
However, an important role is played also by the fact 
that an increase (decrease) of the central density of an 



equilibrium WD leads to a practically exact homologous^ 
contraction (expansion) of the WD structure in the in- 
ner regions {m{Tu) < 1 Mq) in the nonrotating case (this 
can be seen in Fig. Ilip . at least in the range of cen- 
tral densities considered in this paper (as we shall see 
in Section [ill Bl this feature also holds to good accuracy 
in the case of rotating WDs). These aspects, in com- 
bination with the homologous nature of WD inner-core 
collapse, make the size and dynamics of the inner core in 
the bounce phase practically independent of the central 
density of t he initial equilibrium WDs. Early analyti- 
cal work [sl. llOOL rilOl | demonstrated (neglecting thermal 
corrections |113| and rotation) that the mass Mfc of the 
inner core is proportional to in the infall phase dur- 
ing which the fluid pressure is dominated by the contribu- 
tion of degenerate electrons. Around bounce, at densities 
near nuclear matter density, the nuclear component dom- 
inates and the simple dependence does not hold ex- 
actly any longer. As discussed in Sec. Ill B 4} we adopt the 
parametrization Ye{p) as extracted from the simulations 
of Dessart et al. [2^ which predict very efficient electron 
capture, resulting in an average inner-core Ye at bounce 
of ^ 0.18 in the high-T models. This is significantly lower 
than in standard iron core collapse where the inner-core 
Yf, at boun ce is expected to be around ^ 0.25 — 0.30 
[iO) [Hi [13] ■ In our nonrotating AIC models, we find in- 
ner core masses at bounce Mic,b ^ 0.27 Mq (see Fig. [T^ 
which are, as expected, significantly sm aller t han in iron 
core collapse (where Mic,b ~ 0.5Mq [H, [HI, [lil|) . Due 
to their small mass, our AIC inner cores have less ki- 
netic energy at bounce and reach lower densities than 
their iron core counterparts. For example, the nonro- 
tating AIC models exhibit central densities at bounce 
of ~ 2.8 X lO^"' g cm~'^, while in nonrotating iron core 
collapse, maximum densities of > 3 x lO^"' g cm~'^ are 
generally reached at bounce in simulations {e.g., [1^). 
In addition to A/ic,b, the bounce density depends also on 
the stiffness of the nuclear E OS w hose variation we do 
not explore here (see, e.g., [iol . Ill6l |). 

Since the free proton fraction at precollapse and early 
collapse densities grows rather rapidly with te mpe rature 
in the range from ~ 10^ to ~ 10^° K {e.g., jll7l |). the 
efficiency of electron capture is sensitive to the temper- 
ature of WD matter. For example, in the low-T mod- 
els, the value of Yp drops to ~ 0.32 when the density 



^ For a d iscussion of homology in the stellar structure context, 
see [iTl . 

* We define the inner core as the region which is in sonic contact 
at the time of bounce, i.e., 

Afi,,b = [ pWdV , (17) 

J\vr\<c^ 

where W is the Lorentz factor and dV is the invariant 3-volume 
element. The bounce time is defined as the time when the radial 
velocity of the outer edge of the inner core becomes positive. 
Note that such a measure of the inner core is strictly vaUd only 
at the time of bounce. 
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reaches lO^^gcm"'^, while in the high-T models, we ob- 
tain Ye ~ 0.3 at that time. Due to this dependence of Ye 
on T, the inner core masses of low-T models are larger by 
~ 10% compared to those of high-T models. Moreover, 
since the electron degeneracy pressure is proportional to 
(Yep)^^^ [1], the collapse times of the low-T models are 
longer by ~ 5%. We find similar systematics in test calcu- 
lations in which we modify the Ye{p) trajectories of low-T 
models to yield larger Ye at bounce (see Sec. Ill C 3)) . An 
increase of the inner-core Ye by 10% (20%) leads to an 
increase of Mi^.h by 11% (~ 25%). 

It is important to note that the nonrotating AIC mod- 
els discussed above as well as all of the other models 
considered in this study, experience prompt hydrody- 
namic explosions. The bounce shock, once formed, does 
slow down, but never stalls and steadily propagates out- 
wards. While the shock propagation is insensitive to the 
initial WD temperature profile, it shows significant de- 
pendence on the initial WD central density: Owing to 
the greater initial compactness and the steeper density 
gradient of the higher-density models, the shock propa- 
gation in those models is faster and the shock remains 
stronger when it reaches the WD surface. For example, 
in the lowest-density progenitor model AUO, the shock 
reaches the surface within ~ 120 ms after its formation, 
while in the highest density model DUO, it needs only 
~ 80 ms. We point out that Dessart et al. [H, [2^ and 
previous AIC studies [2^, [23| reported significant shock 
stagnation in the postbounce phase of AIC due to the dis- 
sociation of infalling material and neutrino losses from 
the postshock region. Our present computational ap- 
proach includes dissociation (through the EOS, see Sec. 
Ill A 1[) . but does not account for neutrino losses in the 
postbounce phase. Hence, the "prompt" explosions in 
our models are most likely an artifact of our incomplete 
treatment of the postbounce physics. 



B. Rotating AIC 

The AIC of rotating models proceeds through the 
same stages as AIC without rotation and exhibits sim- 
ilar general features, including the well defined split of 
the WD into an inner core that is in sonic contact and 
collapses quasi- homologously^, and a supersonically in- 
falling outer core. Conservation of angular momentum 
leads to an increase of the angular velocity fl oa zu~'^ 
and of the centrifugal acceleration acont = f^^n? oc w~^. 
The latter has opposite sign to gravitational acceleration, 
hence provides increasing centrifugal support during col- 
lapse, slowing down the contraction and, if sufficiently 
strong, leading to centrifugally-induced core bounce only 



^ The collapse is ^uasi-homologous because in this case the rela- 
tion between the infall velocity depends on both the radial 
coordinate and on the polar coordinate [TsI . 
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FIG. 8: Times to core bounce from the onset of collapse as 
a function of initial central angular velocity flc,i- Shown are 
the results of the high-T model sequence (the low-T mod- 
els exhibit identical systematics). Models denoted by an un- 
filled (filled) circle undergo a pressure-dominated bounce with 
(without) significant prompt postbounce convection. Models 
marked by a cross undergo centrifugal bounce at subnuclear 
densities and models marked with a small (large) symbol are 
of set A (D). The colors correspond to various precollapse ro- 
tational configurations (see the legend in Fig. [9]). Note that 
due to their much higher initial compactness, the high-density 
D models (shown in the inset plot) have much shorter collapse 
times than their lower-density A counterparts. 
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FIG. 9: The maximum density pmax,b at bounce as a function 
of the inner core parameter /?ic,b at bounce for the entire set of 
high-T AIC models. Due to the increasing role of centrifugal 
support, pmax.b dccreascs monotonically with increasing rota- 
tion (see the main text for details). The symbol convention 
for the various sets is explained in the caption of Fig. [8] 
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FIG. 10: The inner core parameter /3ic,b at the time of bounce 
for all low-T AIC models plotted as a function of the prec- 
ollapse central angular velocity fie,!- For models with slow 
to moderately-rapid rotation, /9ic,b increases roughly linearly 
with r2c,i (at fixed pc,i and rotation law). In more rapidly 
rotating models, the growth of /3ic,b saturates at ~ 24.5%, 
and further increase of the progenitor rotation results in a 
decrease of /3ic,b. Since D models experience less spin-up dur- 
ing collapse than A models, an increase of at fixed f2c,i 
and rotation law results in a decrease of /3ic,b. A uniformly 
rotating model with a given fie,! and pc,i reaches smaller Ac,b 
than the differentially rotating model with the same f2c,i and 
Pc,i. The symbol convention for the various sets is explained 
in the caption of Fig. [S] 



sligh tly above nuclear density or even at subnuclcar den- 
sity [l08l . flT8l |. 

Just as in the case of nonrotating AIC, models of set A 
collapse more slowly than D models because of the depen- 
dence of the collapse times on the initial central densities. 
However, due to centrifugal support, the collapse times 
grow with increasing precollapse rotation. This is visual- 
ized in Fig.[5]in which we plot the time to core bounce as 
a function of the initial central angular velocity Sl^i. The 
maximum angular velocity of uniformly rotating models 
is limited by the WD surface mass-shedding limit and is 
- 3.5 rad s"^ (- 9.5rads"^) in model AU5 (DU7). The 
effect of rotation on the collapse time of uniformly rotat- 
ing models (dark-green symbols in Fig. is small and 
the time to core bounce increases by ^ 5% from zero 
to maximum precollapse rotation in model set A. The 
more compact D models collapse much faster than their 
lower initial density A counterparts and, in addition, ex- 
perience a smaller spin-up of their more compact inner 
cores. Hence, uniformly rotating D models are less af- 
fected by rotation and their collapse times vary by only 
~ 0.8% from zero to maximum rotation. 

As mentioned in Sec. Ill B 21 WD models that rotate 
differentially according to the rotation law of Yoon & 
Langer 0, [23j have an angular velocity that increases 
from its central value ric,i with tu up to a maximum 



^^max,i at the cylindrical radius tUp, beyond which 17 de- 
creases to sub-Keplerian values at the surface (see Fig.[3|). 
The rate at which increases in the WD core is con- 
trolled by the shear parameter /sh, which we choose in 
the range from 0.2 to 1. The case /sh = 0.2 corre- 
sponds to a nearly uniformly rotating inner region, while 
/sh = 1 corresponds to strong differential rotation with 
il(tUp)/ilc.i ~ 2 — 3. In mass coordinate, this corre- 
sponds to a ratio ^l{Mic^h)/^c,i of ~ 1.4 — 2.4, where 
A/ic.b — 0.3 Mq is the approximate mass that constitutes 
the inner core at bounce in a nonrotating WD model. 

In contrast to uniformly rotating models, differentially 
rotating WDs are not limited by the mass shedding limit 
at the surface. As a result, fic.i can in principle be in- 
creased up to the point beyond which the precollapse WD 
inner core becomes fully centrifugally supported and does 
not collapse at all. For model set AD, this maximum of 
flcA is ^ 8rads"^ (the low-T model AD13f4, which be- 
comes centrifugally supported already at a central den- 
sity of ^ 7 X 10^°gcm~^) while the more compact DD 
models still collapse rapidly at Sic,! ~ ISrads"^ (model 
DD7) . As shown in Fig. [U the most rapidly rotating AD 
model (AD13f4) reaches core bounce after a time which 
is ^ 55% larger than a nonrotating A model. For the 
most rapidly rotating DD model this difference is only 
-5%. 

In Fig. [S] we plot the maximum density /Omax,b at 
bounce as a function of the inner core parameter /3ic.b 
at bounce. Slowly to moderately-rapidly rotating WDs 
that reach Pi^h ^ 15% are only mildly affected by ro- 
tation and their /Omax.b decrease roughly linearly with 
increasing /3ic,b, but stay close to Pnuc- The effect of ro- 
tation becomes nonlinear in more rapidly rotating WDs. 
Models of our set that reach /3ic,b ^ 18% {i.e., AD mod- 
els with n^i ^ 5rads~^) undergo core bounce induced 
partly or completely centrifugally at subnuclear densities. 

As shown in Fig. 1101 Ac.b is a monotonic function 
of Qc.i, but is very sensitive to both the rotation law 
and the initial WD compactness. Our most rapidly uni- 
formly rotating models AU5 and DU7 (both near the 
mass-shedding limit) reach /3ic,b of ~ 10.5% and ~ 3.6%, 
respectively. Hence, uniformly rotating WDs always un- 
dergo core bounce due to the stiffening of the nuclear 
EOS and with little influence of rotation on the dynam- 
ics. 

In models where centrifugal effects remain subdom- 
inant during collapse, /?ic,b grows practically linearly 
with ric,i. This relationship flattens off for models that 
become partially or completely centrifugally supported 
near bounce. Pic.h grows with increasing rotation up to 
~ 24.5% (model AD13f4), beyond which any further in- 
crease in precollapse rotation leads to a decrease of Pic.h, 
since the inner core becomes fully centrifugally supported 
before reaching high compactness and spin-up. In other 
words, there exists a "centrifugal limit" beyond which 
centrifugal forces dominate and, as a result, increasing 
precollapse rotation leads to a decreasing /?ic,b at core 
bounce. This result is analogous to what previous stud- 
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FIG. 11: Mass Mo(rc) in units of solar masses of the WD inner 
region plotted as a function of the equatorial radial coordinate 
re for a number of AIC models with slow and rapid rotation as 
well as high and small central densities. The inset plot shows 
the same on a linear radial scale. The mass distribution of 
the inner Mo(ro) < 0.5Mq region is largely independent of the 
rotational configuration, while an increase (decrease) of the 
central density leads to homologous contraction (expansion) 
of the inner regions. See Sec. Ill CI for details of the model 
setups. 
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FIG. 12: Angular velocity profiles in the equatorial plane at 
the time of core bounce and time-evolution of the central den- 
sity pc (in the inset) for models AU2 and DU7. The initial 
angular velocity of model DU7 is larger by a factor of ^ 5.3 
than that of model AU2, but the latter experiences a ~ 5.3 
greater spin-up during collapse. As a result, these models 
produce inner cores with almost identical rotational config- 
urations and similar masses in the bounce phase. This is 
refiected in an identical evolution of the central densities at 
bounce. 



ies [i^ llOTf found in the rotating core collapse of massive 
stars and has consequences for the appearance of nonax- 
isymmetric rotational instabilities in PNSs. This will be 
discussed in more detail in Sec. |Vl 

The influence of the precollapse compactness on the 
dynamics of rotating AIC can also easily be appreciated 
from Fig. (TUl The higher-density, more-compact WDs of 
set D spin up much less than their A counterparts since 
their inner cores are already very compact at the onset 
of collapse. Hence, a higher-density WD that reaches a 
given value of /3ic,b must have started out with a larger 
J7c,i than a lower-density WD reaching the same /3ic,b- 
For the particular choice of initial central densities rep- 
resented by D and A models and in the case of uniform 
or near-uniform rotation, the ratio between the f2c,i of a 
D and A model required to reach the same /3ic,b is ~ 5.3. 
This factor can be understood by considering Fig. [11] in 
which we plot the enclosed mass as a function of equa- 
torial radial coordinate Vc of selected A and D initial 
WD configurations with slow and rapid rotation. The 
important thing to notice is that the WD core structure 
{M < O.SMq) is insensitive to the rotational configura- 
tion and obeys a homology relation. Stated differently, 
for a model of set D, a homologous expansion in the radial 
direction by a factor of ^ 2.3 yields an object whose inner 
part is very similar to a lower-density A model. In turn, 
the collapse of A models corresponds to a ~ 2.3 times 
greater contraction of the WD core compared to their D 
model counterparts and a spin up that is greater by a 
factor of ^ (2.3)^ ~ 5.3. This explains the strong depen- 



dence of the inner-core angular velocity and Ac,b on the 
initial central density observed in Fig. 1101 Furthermore, 
it suggests that one can find A— D model pairs that differ 
greatly in their precollapse angular velocities, but yield 
the same rotational configuration at bounce. An exam- 
ple for this is shown in Fig. [T^] in which we plot for the 
uniformly rotating model pair AU2-DU7 the equatorial 
angular velocity profile at the time of bounce as well as 
the evolution of the central density around the time of 
bounce. AU2 and DU7 have practically identical angu- 
lar velocity profiles and their core structure, core mass 
and /3ic,b agree very closely. As can be seen in the inset 
plot of Fig. [T^l this results in nearly identical pc time 
evolutions around bounce and demonstrates that WDs 
with quite different precollapse structure and rotational 
setup can produce identical bounce and postbounce dy- 
namics. This can also occur for pairs of differentially ro- 
tating models and is an important aspect to keep in mind 
when interpreting the GW signal from AIC discussed in 
SeclIVl 

Figure [T3] shows the mapping between /3ic.b and the 
inner core mass Mic.b at bounce for all high-T models. 
Rapid (differential) rotation not only increases the equi- 
librium mass of WDs (sec Table |l|, but rotational sup- 
port also increases the extent of the region in sonic con- 
tact during collapse. Hence, it may be expected that 
-^^icb grows with increasing rotation. However, for WDs 
below /Jic.b S 13%, Mic.b is essentially unaffected by ro- 
tation and stays within 0.02 Mq of the nonrotating value 
of 0.28 A/q. Only when the effects of rotation become 
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from to 5 radms^^ (~ 3.3 rad ms^^), while their 
pole to equator axis ratios vary from 1 to ^ 0.4 (~ 0.6)^. 
Some of the rapidly rotating WDs produce PNSs with a 
slightly off-center maximum in density, though the den- 
sity distribution of the inner regions does not exhibit a 
pronounced toroidal geometry. The clearest deviation 
from a centrally-peaked density distribution is produced 
in the case of model ADIO, which reaches /?ic,b — 21.3 % 
(/3ic,b ^ '24:A %) in its high-T (low-T) variant. In this 
model, the point of highest density after bounce is located 
at r ~ 0.94 km, but the maximum value is larger than the 
central density value by only ^ 0.3%. For models with 
less rapid rotation, the off-center maximum is much less 
pronounced, and completely disappears for /3ic b below 
-20%. 



FIG. 13: Mass Aficb of the inner core at bounce for all high-T 
models versus the parameter /3ic,b of the inner core at bounce. 
The systematics of Mic,b with Ac,b are identical for the set 
of low-r models, but their Mic,b are generally ~ 10% larger. 
The symbol convention for the various sets is explained in the 
caption of Fig. [S] 



strong at /?ic,b 13 — 18% does Mjc.b increase roughly 
linearly with /3ic,b. WDs that undergo centrifugal bounce 
have /9ic,b ^ 18% and correspondingly large inner cores 
that are more massive than — 0.5 Mq. Such values of 
Mic,b are accessible only to differentially rotating WDs. 

Also for rotating models, the dependence of the AIC 
dynamics on the initial temperature of AIC progenitor 
WDs is simple and straightforwardly understood from 
the nonrotating results discussed in Sec. IIII Al These 
show that the low-T models yield inner cores that are 
~ 10% larger in mass than in their twice-as-hot high-T 
counterparts. Due to their larger mass, the inner cores 
of collapsing low-T AIC progenitors also contain a larger 
amount of angular momentum. At fixed rotation law and 
f2c,i, they reach values of /3ic,b that are larger by up to 
~ 5% (in absolute value!). Hence, lower-T WDs become 
affected by centrifugal support, bounce centrifugally, and 
reach the centrifugal limit at lower than their higher- 
T counterparts. Along the same lines behave test calcu- 
lations in which we impose increased inner-core values of 
Ye (see Sees. Ill A 2] and IIIC3)) . The increased Ye leads 
to more massive and more extend inner cores which, in 
turn, are more likely to experience a centrifugal support. 

To conclude our discussion of rotating AIC, we sum- 
marize for the reader that the PNSs born from the set 
of differentially (uniformly) rotating AIC models consid- 
ered here have average angular velocities^ in the range 



C. Shock Propagation and the Formation of 
Quasi-Keplerian Disks 

As pointed out earlier (see Sec. IIII Ap . all AIC mod- 
els considered in this study undergo weak hydrodynamic 
explosions. This is an artifact of our approach that ne- 
glects postbounce neutrino emission, but is unlikely to 
strongly affect the results presented in this section, since 
in the MGFLD simulations of [l^, the shock stalls only 
for a very short period and a weak explosion is quickly 
initiated by neutrino heating. 

In moderately-rapidly and rapidly rotating AIC (with 
Ac,b ^ 5%), the shock propagation is significantly af- 
fected by centrifugal effects. The material near the equa- 
torial plane of rotating WDs experiences considerable 
centrifugal support, and its collapse dynamics is slowed 
down. As a consequence, the bounce is less violent and 
the bounce shock starts out weaker near the equato- 
rial plane than along the poles. Centrifugal support of 
low-latitude material also leads to reduced postbounce 
mass accretion rates near the equatorial plane, facilitat- 
ing steady propagation of the shock at low latitudes. In 
the polar direction, where centrifugal support is absent, 
the shock propagates even faster due to the steeper den- 
sity gradient and smaller polar radius of the WD. This 
quickly leads to a prolate deformation of the shock front 
in all rotating models and the shock hits the polar WD 
surface much before it breaks out of the equatorial enve- 
lope. This is shown in Fig. 1141 where we plot the equa- 
torial and polar profiles of the radial velocity and spe- 
cific entropy per baryon for model DD7 at various post- 
bounce times. Due to the prolateness of the shock front, 
it breaks out of the polar surface ^130 ms before reach- 
ing the WD's equatorial surface. Moreover, due to the 



The average angular velocity Q of the differentially rotating 
models considered here is computed using the approximation 
f2 = Jic/hc} where Jit, is the inner core angular momentum and 
lie is the (Newtonian) inner core moment of inertia. 



^ The PNS formed in the AIC of WDs is surrounded by hot low 
density material in the early postbounce phase, making it hard to 
define the boundary of the PNS unambiguously. For the present 
rough estimate of the axis ratio, we assume a density threshold 
of lO'^^ gcm^^ to mark the boundary of the PNS. 
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anisotropy of the density gradient and the initial shock 
strength, the specific entropy of the shock-heated ma- 
terial is larger by a factor of 2 — 3 along the polar 
direction. 

The asphericity of the shock front and the anisotropy 
of the shock strength become more pronounced in AIC 
with increasing rotation [H^, IZ^I • As pointed out in Sec- 
tion IIII Al due to the their greater initial compactness 
and thus steeper density gradients, the shock propagates 
faster in D models: In model DDI, for example, the 
shock reaches the surface in the equatorial plane within 

88 ms, while for model ADl, the corresponding time 
is ^ 143 ms. 

Rapid rotation and, in particular, rapid differential ro- 
tation, increases the maximum allowable WD mass. The 
most rapidly uniformly rotating WDs in our model set 
{i.e., models DU7 and AU5) have an equilibrium mass 
of ^ 1.46M0, which is only slightly above Mch in the 
nonrotating limit. Our most rapidly differentially rotat- 
ing WDs (models AD13f2 and DD7), on the other hand, 
reach equilibrium masses of up to ~ 2 AIq. Much of the 
rotationally supported material is situated at low lati- 
tudes in the outer WD core, falls in only slightly during 



collapse, and forms a quasi-Keplerian disk-like structure. 
The equatorial bounce shock is not sufficiently strong to 
eject much of the disk material and "wraps" around the 
disk structure, producing only a small outflow of outer 
disk material at Vr ^ 0.025c. This is in agreement with 
Dcssart ct al. [l^, who first pointed out that rapidly ro- 
tating AIC produces PNSs surrounded by massive quasi- 
Keplerian disk-like structures in the early postbounce 
phase. As recently investigated by Metzger et al. 
(but not simulated here), the hot disk will experience 
neutrino-cooling on a timescale of ~ 0.1 s, dri ving the 
disk composition neutron-rich to reach Yg '^-^ 0.1 |28l. [isj. 
depleting the pressure support and leading to limited 
contraction of the inner parts of the disk. The outer and 
higher-latitude regions ex pan d with a neutrino-driven 
wind [11] . As discussed by [48j , subsequent irradiation of 
the disk by neutrinos from the PNS increases its proton- 
to-neutron ratio, and Ye may reach values as high as 
~ 0.5 by the time the weak interactions in the disk freeze 
out. The disk becomes radiatively inefficient, a-particles 
begin to recombine and a powerful disk wind develops, 
blowing off most of the disk's remaining material. Met- 
zger et al. [13 argue that, depending on disk mass, the 
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TABLE III: Summary of properties of the quasi-Keplerian 
disks formed in the set of AIC models AD, AU, DD and DU. 
ii/disk is the thickness and Re is the equatorial radius of the 
disk, while Maisk is its mass. These quantities are computed 
at the time when the shock reaches the WD surface in the 
equatorial plane. The disk parameters do not vary signifi- 
cantly between the two choices of WD temperature considered 
in this study. 
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Mdisk ^ O.OSA/q^ and is produced in model AU5 which 
rotates near the mass-shedding limit. Since the angular 
velocity of the outer (tu > vup) core of differentially ro- 
tating models is set to reach nearly-Keplerian values (cf. 
Eq. most of the outer WD envelope has substan- 

tial centrifugal support and thus the differentially rotat- 
ing models yield significantly larger 7\/disk. For example, 
model AD4 which has /3ic.b and total angular momentum 
comparable to model AU5 yields a disk mass of ~ 0.2Mq. 

The total mass of the disk and its equatorial radius (the 
disk thickness iJdisk) grow with increasing rotation (see 
Tab. imp . Slowly rotating models such as ADl have little 
centrifugally supported material and acquire spheroidal 
shape soon after bounce, resulting in a disk mass as small 
as ~ O.O3M0. More rapidly rotating models such as 
ADIO produce significantly more strongly flattened disks, 
with i?o ~ 2700 km, i/disk ^ 800 km and a disk mass 
■^disk of ~ 0.8A'/q. Due to the greater initial compact- 
ness of the higher-density D models, their disks are less 
massive and have smaller equatorial radii when compared 
to A models. Hence, when considering two WDs of set A 
and D with the same total angular momentum, the mass 
and equatorial extent of the disk in the D model will be 
smaller by a factor of 1.5 — 2. 

These results indicate that massive disks of Mdisk ^ 
O.IM0 are unlikely to be compatible with the assump- 
tion of uni formly rotating accreting WDs argued for 
by [l^, Il0l| . In order to produce disks of appreciable 
mass and significant ^^Ni-outflows in AIC, the progeni- 
tor must either be an accreting WD obeying a differential 
rotation law similar to that proposed by 0, , or may 
be the remnant of a binary- WD merger event. However, 
for the latter, the differential rotation law is unknown 
and may be very different from what we consider here 
(for a dis cussion of binary- WD merger simulations, see, 
6.17., [ll9l | and references therein). 



IV. GRAVITATIONAL WAVE EMISSION 



outflows synthesize of the order of lO"'^ — 10"^ Mq of 
^^Ni, but very small amounts of intermediate-mass iso- 
topes, making such AIC explosions spectroscopically dis- 
tinct from ^^Ni outflows in standard core-collapse and 
thermonuclear SNe. 

Our results, summarized in Tab. IIIIl show that the 
masses and the geometry of the disks produced in AIC 
are sensitive to the angular momentum distribution in 
the precollapse WDs. In models with uniform rotation 
below the mass-shedding limit, only a very small amount 
of low-latitude material rotates at near-Keplerian angu- 
lar velocities. Therefore, most of the outer core mate- 
rial of such models undergoes significant infall, so that 
uniformly rotating WDs will generally produce small 
disks. The largest disk mass for uniform rotation is 



In Fig. [ini we present the time evolution of the GW 
strain h at an assumed source distance of 10 kpc for a 
representative set of AIC models evolved with the Ye ip) 
parametrization obtained from (28l |. The GW signals of 
all models have the same overall morphology. This gen- 
eral AIC GW signal shape bears strong resemblance to 
GW signals that have been classified as "type HI" in the 
past [m, [tI, [t^ , but also has features in common with 



* We point out that because the disks do not settle down to exact 
equilibrium right after bounce or not even after shock passage, it 
is hard to introduce an unambiguous definition of the disk mass. 
In the present study, we define the disk as the structure that sur- 
rounds the PNS at ro > 20 km with densities below 10^^ g cm" ^ 
and angular velocity Q > 0.58 Hk. The latter condition ensures 
that the disk cannot contract by more than a factor of ~ 3 as a 
result of cooling. 
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FIG. 15: Evolution ol the dimensionless GW strain h (in 
units of 10~^^ at a source distance of 10 kpc) as a function 
of postbounce time for representative models with different 
precollapse rotation profiles, central densities and tempera- 
tures (low-T models with solid black lines and high-T with 
dashed red lines). Models with slow and (almost) uniform pre- 
collapse rotation {e.g., AD 114 or DU3) develop considerable 
prompt postbounce convection visible as a dominating lower- 
frequency contribution in the waveform. Centrifugal effects 
damp this prompt convection and the waveforms of models 
with moderately-rapid rotation {e.g., ADl, AU4, DD7 and 
ADS) and of rapidly rotating models {e.g., AD9 or ADIO) 
exhibit no such contribution to the signal. 



the GW signals predicted for rotating iron core collapse 
("type r,m- 

The GW strain h in our AIC models is positive in 
the infall phase and increases monotonically with time, 
reaching its peak value in the plunge phase, just ~ 0.1 ms 
before bounce. Then, h rapidly decreases, reaching a 
negative peak value within ~ 1 — 2 ms. While the first 
positive peak is produced by rapid infall of the inner core, 
the first negative peak is caused by the reversal of the 
infall velocities at bounce. Following the large negative 
peak, h oscillates with smaller amplitude with a damping 
time of ~ 10 ms, reflecting the hydrodynamical ringdown 
oscillations of the PNS. Although all AIC models of our 



baseline set produce Type III signals, we can introduce 
three subtypes whose individual occurrence depends on 
the parameter /3ic,b of the inner core at bounce: 

Type Ilia. In slowly rotating WDs (that reach 
Ac.b ^ 0.7%), strong prompt convective overturn de- 
velops in the early postbounce phase, adding a lower- 
frequency contribution to the regular ringdown signal 
{e.g., models ADlf4, DU3). The largest-amplitude part 
of this GW signal type comes from the prompt con- 
vection. Nevertheless, the GW signal produced by the 
bouncing centrifugally-deformed inner core is still dis- 
cernible, with the first positive peak being generally 
larger by a factor of ~ 2 than the first negative peak. 
Subsequent ringdown peaks are smaller by a factor of 
> 3. Wc point out that the observed prompt postbounce 
convection is most likely overestimated in our approach, 
since we do not take into account neutrino losses and 
energy deposition by neutrinos in the immediate post- 
shock region, whose effect will quickly smooth out the 
negative entropy gradient left behind by the shock and 
thus significantly damp this early convective instability 
in full postbounce radi ation- hydrodynamics calculations 
(see, e.g., ^MB)- 

Type Illb. In moderately rapidly rotating WDs that 
reach 0.7 < /3i 



ic,b 



< 



18 % and still experience a pressure- 
dominated bounce, convection is effectively suppressed 
due to a sufficiently lar ge po sitive specific angular mo- 
mentum gradient {e.g., [12lj). Hence, there is no no- 
ticeable convective contribution to the postbounce GW 
signal (see, e.g., models ADl, AU4, DD7, ADS). For 
this signal subtype, the peak GW strain |ft.|max is asso- 
ciated with the first positive peak while relative values 
of the amplitudes of the first several peaks are similar to 
type Ilia. 

Type IIIc. If rotation is sufficiently rapid and leading 
to /3ic.b 18 %, the core bounces at subnuclear densities 
due to strong centrifugal support. This is reflected in the 
GW signal by an overall lower-frequency emission and 
a significant widening of the bounce peak of the wave- 
form (see, e.g., models AD9, ADIO). In some models of 
this subtype, the negative peak can be comparable to or 
slightly exceed that of the first positive peak in the wave- 
form. This reflects the fact that the plunge acceleration is 
apparently reduced more significantly by rotation than is 
the re-expansion acceleration at core bounce. The post- 
bounce ringdown peaks in all type IIIc models are smaller 
by a factor of > 2 compared to the bounce signal. As 
pointed out in Sec. IIII Bl uniformly rotating models do 
not rotate sufficiently rapidly to experience centrifugal 
bounce. Hence, they do not produce a type IIIc signal. 

The AIC GW signal morphology is affected only 
slightly by variations in WD temperature and their re- 
sulting changes in the inner core that are on the few- 
percent level for the range of precollapse temperatures 
considered here. In test calculations with more substan- 
tially increased inner-core values of Ye (sec Sees. Ill A 21 
and III C 3p and, in turn, significantly larger values of 
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FIG. 16: Peak value |/i|max of the GW amplitude at a source 
distance of 10 kpc distance for the burst signal of all models 
versus the parameter /3ic,b of the inner core at the time of 
bounce. At slow to moderately rapid rotation, |/!.|max scales 
almost linearly with /3ic,b (marked by the dotted straight line) , 
while for /?b > 16% centrifugal effects reduce |/i|max. The 
symbol convention for the various sets is explained in the 
caption of Fig. [S] 



.^^ic,bj wc find signals that arc intermediate between 
type III and type I. 

Key quantitative results from our model simulations 
are summarized in Tables . HIl and HVl The wave form data 
for all models are available for download from Il22l. 



A. Peak Gravitational Wave Amplitude 

Across our entire model set, the peak GW amplitude 
l^lmax covers a range of almost two orders of magnitude, 



from ~ 10 



to 



10" 



(at distance to the source of 



10 kpc; see Tab. HI)) . |/i|max depends on various parame- 
ters and it is difficult to provide a simple description of 
its systematics that encompasses all cases. In order to 
gain insight into how |/i|max depends on f2c,i, on differ- 
ential rotation, on the initial pc, on the precollapse WD 
temperature, and on the degree of deleptonization in col- 
lapse, we describe below the effects of variations in one 
of these parameters while holding all others fixed. 

(i) In a sequence of precollapse WDs with fixed dif- 
ferential rotation, pc, and Tq, the peak GW amplitude 
l^lmax increases steeply with in slowly rotating mod- 
els that do not come close to being centrifugally sup- 
ported. When centrifugal effects become dynamically 
important, |ft.|max saturates at ~ 7 x 10~^^ (at 10 kpc) 
and then decreases with increasing r2c,i- This reflects the 
fact that such rapidly spinning inner cores produced by 
AIC cannot reach high densities and high compactness 
and that the slowed-down collapse decreases the decel- 
eration at bounce, thus reducing |ft.|max and pushing the 
GW emission to lower frequencies. 




FIG. 17: Spectral energy density of the GW signal for repre- 
sentative AIC models AD0f3 (top panel), AD4 (center panel), 
and ADIO (bottom panel), /max is the peak frequency of the 
GWs emitted at core bounce. 



(ii) In a sequence of precollapse WDs with fixed fJc,i, 
To, and pc, an increase in the degree of WD differen- 
tial rotation leads to an increase in the amount of angu- 
lar momentum present in the WD inner core at bounce. 
This translates into an increase of |/i|max in models that 
do not become centrifugally supported and experience a 
pressure-dominated bounce. The transition to centrifu- 
gal bounce is now reached at lower values of Qc,i (see 
Sec. nil Bp . so that the centrifugal saturation of |/i|max 
described above in (i) is reached at much smaller values 

of r^ci- 

(Hi) In a sequence of precollapse WDs with fixed ^c,i, 
Tq, and differential rotation and varying pc, models with 
lower (higher) pc yield larger (smaller) values of |/i|max- 
This is because models that are initially less compact spin 
up more during collapse (c/. the discussion in Sec. IIII Bj) . 
However, this systematics holds only as long as the model 
does not become centrifugally supported, which happens 
for lower (higher) pc WDs at smaller (greater) ilc.i- 

(iv) When only the WD temperature is varied, we 
find that for slowly to moderately-rapidly rotating WDs, 
high-T models generally reach smaller |/i|max than their 
low-T counterparts. This is because high-T WDs yield 
smaller inner cores at bounce, which hold less angular 
momentum and, as a consequence, are less centrifugally 
deformed (see Tab. HIl and in Fig. [T5l) . However, this be- 
havior reverses in rapidly rotating WDs for which low-T 
models are more centrifugally affected and, hence, yield 
a smaller |/i|max than their high-T counterparts. 

(v) If the degree of deleptonization is decreased by 
an ad-hoc increase of inner-core Ye (see Sees. Ill A 21 and 
III C 3p and all else is kept fixed, Mic,b increases and for 
slowly to moderately rapidly rotating WDs, |/i|max in- 
creases. As for the low-T case discussed in the above, 
this behavior reverses in rapidly rotating WDs for which 
high- Ye models are more centrifugally affected and yield 
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FIG. 18: Frequency /max at the maximum of the GW spectral 
energy density in pressure-dominated bounce and a subset of 
centrifugal bounce models versus the parameter /3ic,b of the 
inner core at bounce. The symbol convention for the various 
sets is explained in the caption of Fig. [S] 

smaller |/i|max than their lower- l^e counterparts. 

To demonstrate the dependence of |ft.|max on the over- 
all rotation of the inner core at bounce, we plot in Fig. [TBI 
l^lmax as a function of the inner core parameter /3ic,b at 
bounce for our high-T models. |/i|max depends primar- 
ily on /3ic,b and is rather independent of the particular 
precollapse configuration that leads to a given /3ic,b- For 
small Ac,b far away from the centrifugal limit, we find 
l^lmax oc /Sjc^tf, whcrc we have obtained the exponent by 
a power-law fit of high-T models with 1% < /Jjc^b ^ 13%. 
This finding is in qualitative agreement with what [49| 
saw for iron core collapse. The overall maximum of |ft.|max 
is reached in WDs that yield /3ic,b ^ 16%, beyond which 
l^^lmax decreases with increasing /3ic,b. 

B. Gravitational- Wave Energy Spectrum 

The total energy emitted in GWs is in the range of 
~ 10-1'^ Mqc^ < £:gw < 2 X lO-^Mgc^ in the entire set 
of models considered in this article. In Fig. [T71 we plot 
the GW spectral energy density dEa^/df of the three 
models ADlfS, AD4 and ADIO as representative exam- 
ples of the three signal subtypes Ilia - IIIc. The top 
panel shows model ADlfS as a representative pressure- 
dominated bounce model with prompt convection. In 
such models, there is a strong structured, but broad con- 
tribution to the spectrum at low frequencies. The inte- 
gral of such a contribution (which is present in all models 
with slow rotation) can exceed that from core bounce in 
these models. This is not the case in model ADlfS, whose 
GW burst from bounce is the one leading to the peak at 
/„,ax = 720 Hz. 

The central panel of Fig. [17] depicts dEQy^/df of model 
AD4 as a representative pressure-dominated bounce 
model in which no significant postbounce convection oc- 
curs. The spectrum of this model exhibits a distinct 
and narrow high-frequency peak at /max ~ 805 Hz. Fi- 
nally, the bottom panel of Fig. [T7| refers to model AlO 



that experiences centrifugal bounce. In this model, the 
dynamics is dominated by centrifugal effects, leading to 
low frequency emission and /max = 310 Hz, but higher- 
frequency components are still discernible and arc most 
likely related to prolonged higher-frequency GW emission 
from the PNS ringdown. 

In Fig.[T8l we plot the peak frequencies /max of the GW 
energy spectrum as a function of the inner core parameter 
/?ic,b for high-T AIC models (the low-T and higher- Y; 
models show the same overall systematics). In models 
that undergo pressure-dominated bounce, /max increases 
nearly linearly with /3ic,b in the region Pic,h ^ 10 %, while 
at /3ic,b in the range of 10 % < /3ic,b ^ 20 %, the growth 
of /max saturates at ~ 800 Hz and /max does not change 
significantly with further increase of rotation. For very 
rapid rotation (/3ic,b ^ 20 %), /max decreases steeply with 
/?ic,b, reaching a value of ~ 400 Hz at (3ic,h — 23% (not 
shown in the figure, see Table |TV|. 

While it is straightforward to understand the system- 
atics of /max at high /3ic,b where centrifugal effects slow 
down collapse and thus naturally push the GW emission 
to low frequencies, the increase of /max with rotation 
at low to intermediate /3ic,b is less intuitive. If one as- 
sumes that the dominant GW emission at core bounce 
in all models is due to the quadrupole component of 
the fundamental quasi-radial mode of the inner core, one 
would expect a monotonic decrease of /max with increas- 
ing rotation and, hence, decreasing mean core density 
(see, e.g., [i^). A possible explanation for the increase 
of /max at slow to modcrately-rapid rotation is that the 
primary GW emission in these models is due to the fun- 
damental quadrupole ^/-modc, whose frequency may in- 
crease with rotation^ This has been demonstrated by 

Dimmelmeicr et al. |l23l | who studied oscillation modes 
of sequences of 7 = 2 polytropes. To confirm this inter- 
pretation, a nd fo llowing the technique of mode-recycling 
outlined in |l23l | , we perturb a subset of our postbounce 
cores with the eigenfunction of the ^/-mode of a New- 
tonian nonrotating neutron star. As expected, we find 
that the resulting dynamics of the postbounce core is 
dominated by a single oscillation mode with a frequency 
that matches within < 10% the peak frequency /max 
of dEQy//df of the corresponding slowly or moderately 
rapidly rotating AIC model. The interesting details of 
the mode structure of the inner cores of AIC and iron 
core collapse will receive further scrutiny in a subsequent 
publication. 

Finally, in Fig. I19i we provide time-frequency analy- 
ses of the GW signals of the same representative models 
shown in Fig. 1171 The analysis is carried out with a 
short-time Fourier transform employing a Gaussian win- 
dow with a width of 2 ms and a sampling interval of 
0.2 ms. In all three cases, the core bounce is clearly vis- 
ible and marked by a broadband increase of the emit- 
ted energy. The slowly-rotating model ADlf3 emits its 
strongest burst at 600 - 800 Hz (/max = 720 Hz) and 
subsequently exhibits broadband emission with signifi- 
cant power at lower frequencies due to prompt convec- 
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FIG. 19: Time-frequency colormaps of the GW signals of models ADlfS (type Ilia), AD4 (type Illb), and ADIO (type IIIc, see 
also Fig. I17|l . Plotted is the "instantaneous" spectral GW energy density dEow/df in a 2-ms Gaussian window as a function 
of postbounce time. Note the prebounce low-frequency contribution in the moderately-rapidly rotating models (model AD4, 
center panel) and rapidly rotating models (model ADIO, right panel). The range of the colormap of the left panel (model 
ADlfS) is smaller by one dex than those of the other panels. 



tion. Model AD4 is more rapidly rotating and shows 
significant pre-bounce low-frequency emission due to its 
increased rotational deformation. At bounce, a strong 
burst, again with power at all frequencies, but primar- 
ily at frequencies about its /max = 805 Hz, is emitted. 
Much of the postbounce i^cw is emitted through ring- 
down oscillations at /max that may be related to the ^/ 
mode of this model's PNS. Finally, in the rapidly rotat- 
ing and centrifugally bouncing model ADIO, we observe 
again low-frequency emission before bounce, but only a 
small increase of the primary emission band at bounce to 
~ 200 — 400 Hz. Nevertheless, there is still an appreciable 
energy emitted from higher-frequency components of the 
dynamics at bounce and postbounce times. 



Comparing GW signals from AIC 
and Iron Core Collapse 



Recent studies [49|, [6^ [631, l78[ have shown that the col- 
lapse of rotating iron cores (ICC) produces GW signals 
of uniform morphology (so-called "type I" signals, see, 
e.g., [t^) that generically show one pronounced spike as- 
sociated with core bounce with a subsequent ringdown 
and are similar to the type HI signals found here for 
AIC. As in the AIC case, the GW signal of ICC has sub- 
types for slow, moderately rapid, and very rapid rota- 
tion. For comparing AIC and ICC GW signals, we chose 
three representative AIC models and then pick three 
ICC models with similar /3ic,b from the study of Dim- 
melmeier et al. [i^ whose waveforms are freely available 
from |124| . This should ensure that we compare collapse 
models that arc similarly affected by centrifugal effects 
for a one-to-one comparison. However, one should keep 



in mind that the inner core masses Mjc.b at bounce of 
ICC models are generally larger by ^ 0.2 — 0.3 Af0 than 
in our AIC models (see Sec. Illli^ . 

In Fig.[2niwc present this comparison and plot the GW 
signals of the high-T AIC models DU2 (slow rotation, 
type Ilia), AU5 (moderately-rapid rotation, type Illb), 
and ADIO (very rapid rotation, type IIIc). In the 
same order, we superpose the GW signals of the Dim- 
mclmeier et al. [H] models s20A1005, s20A3O07 and 
s20A3O15. These models started with the precollapse 
iron core of a 20 M© star and were run with the same 
code, EOS, and deleptonization algorithm as our AIC 
models, though with different, ICC specific Ye{p) trajec- 
tories. 

The left panel of Fig. [201 compares the slowly rotat- 
ing models DU2 and s20AlO05 that undergo pressure 
dominated bounce and exhibit strong postbounce con- 
vection. As pointed out before, the latter is most likely 
overestimated in our current approach as well as in Dim- 
melmeier et al.'s. Note that the width of the waveform 
peaks associated with core bounce is very similar, indicat- 
ing very similar emission frequencies. Model s20A1005 
exhibits a significantly larger signal amplitude at bounce. 
This is due to s20AlO05's larger Mic.h but also to the fact 
that its Ac.b is ~ 0.7% compared to the ~ 0.4% of DU2 (a 
closer match was not available from |124| |'). The prompt 
convection in model s20A1005 is more vigorous and gen- 
erates a larger- amplitude GW signal than in model DU2. 
This is due to the much steeper density gradient in the 
WD core that allows the AIC shock to remain stronger 
out to larger radii. Hence, it leaves behind a shallower 
negative entropy gradient, leading to weaker convection 
and postbounce GW emission. 

In the central panel of Fig. [201 we compare two 
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FIG. 20: Evolution of GW signals in the high-T AIC models DU2, AU5, ADIO (black solid lines) and three massive star iron 
core-collapse models s20A1005, s20A3O07 and s20A3O15 (red dashed lines) from the model set of Dimmelmeier et al. (igt - 
The AIC model DU2 and the iron core collapse model s20A1005 undergo pressure-dominated bounce with significant prompt 
convection. Models AU5 and s20A3O7 experience pressure-dominated bounce without significant convection, and models 
ADIO and s20A3O15 undergo centrifugal bounce. The inner cores of models DU2, AU5 and ADIO (s20A1005, s20A3O07 and 
S20A3O15) reach values of Ac.b of about 0.4%, 10.2% and 21.3% (0.7%, 10.1%, and 21.6%). Times are given relative to the 
time of core bounce which we mark with a vertical line. 



moderately-rapidly rotating models with nearly identical 
/3ic of ~ 10%. Both models show a pre-bounce rise due to 
the inner core's accelerated collapse in the plunge phase. 
The AIC inner core, owing to its lower Ye and weaker 
pressure support, experiences greater acceleration and 
emits a higher-amplitude signal than its ICC counterpart 
in this phase. At bounce, the stilf nuclear EOS deceler- 
ates the inner core, leading to the large negative peak in 
the GW signal. Because of the more massive inner core in 
ICC and since the EOS governing the dynamics is identi- 
cal in both models, the magnitude of this peak is greater 
in the ICC model. Following bounce, the ICC model's 
GW signal exhibits a large positive peak of comparable 
or larger amplitude than the pre-bounce maximum. This 
peak is due to the re-contraction of the ICC inner core 
after the first strong expansion after bounce. With in- 
creasing rotation, this re-contraction and the associated 
feature in the waveform become less pronounced. On the 
other hand, due to its smaller inertia, the AIC inner core 
does not significantly overshoot its new postbounce equi- 
librium during the postbounce expansion. Hence, there 
is no appreciable postbounce re-contraction and no such 
large positive postbounce peak in the waveform. 

Example waveforms of AIC and ICC models experienc- 
ing core bounce governed by centrifugal forces are shown 
in the right panel of Fig. [501 In this case, the pre-bounce 
plunge dynamics are significantly slowed down by cen- 
trifugal effects and the GW signal evolution is nearly 
identical in AIC and ICC. At bounce, the more mas- 
sive inner core of the ICC model leads to a larger and 
broader negative peak in the waveform and its ringdown 
signal exhibits larger amplitudes than in its AIC coun- 
terpart. 

Finally, we consider AIC models with variations in the 
inner-core Yg due either to different precollapse WD tem- 
peratures or ad- hoc changes of the Y^ (p) parametrization 
(see Sees. Ill A 21 and IIIC 3|) . Lowcr-T WDs yield larger 



inner-core values of Kg and, in turn, larger Mjc.b and 
GW signals that are closer to their iron-core counter- 
parts. The same is true for models in which we impose an 
increased inner-core Y^: AIC models with inner-core Yg 
10% larger than predicted by HI] still show clear type III 
signal morphology while models with 20% larger Ye fall 
in between type III and type I. 

To summarize this comparison: rotating AIC and ro- 
tating ICC lead to qualitatively and quantitatively fairly 
similar GW signals that most likely could not be distin- 
guished by only considering general signal characteristics, 
such as maximum amplitudes, characteristic frequencies 
and durations. A detailed knowledge of the actual wave- 
form would be necessary, but even in this case, a distinc- 
tion between AIC and ICC on the basis of the comparison 
presented here would be difficult. It could only be made 
for moderately-rapidly spinning cores based on the pres- 
ence (ICC, type lb) or absence (AIC, type Illb) of a first 
large positive peak in the waveform, but, again, only if 
AIC inner cores indeed have significantly smaller Y^ than 
their iron-core counterparts. ICC and AIC waveforms of 
types la/IIIa and Ic/IIIc are very similar. Additional 
astrophysical information concerning the distance to the 
source and its orientation as well as knowledge of the 
neutrino and electromagnetic signatures will most likely 
be necessary to distinguish between AIC and ICC. 



D. Detection Prospects for the Gravitational Wave 
Signal from AIC 

In order to assess the detection prospects for the GW 
signal from AIC, we evaluate the characteristic signal fre- 
quency /c and the dimensionless characteristic GW am- 
plitude he- Both quantities are detector-dependent and 
are computed using Eq. and respectively. 

In Fig.inU we show he for all models as a function of /c 
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FIG. 21: Detector-dependent characteristic amplitudes of the 
GW signals of all models at an assumed distance of 10 kpc. 
The symbol convention for the various sets is described in the 
caption of Fig. [S] See the text for a discussion of the numbers 
and arrows. 



for an initial 4-km LIGO detector, assuming a source dis- 
tance of 10 kpc. For comparison with detector sens itivity, 
we include initial LIGO's design h^ms curve [l25j |. The 
distribution of our set of models in this figure obeys sim- 
ple systematics. A number of very slowly rotating mod- 
els that undergo pressure-dominated bounce with prompt 
convection (type Ilia) form a cluster in frequency in one 
region (near arrow 1). These models have the overall low- 
est values of he and exhibit low values of fc in the range 
of 130 — 350 Hz. Both fc and he grow with increasing 
rotation (along arrow 1). For the pressure-dominated 
bounce models without significant prompt convection 
(type Illb), he grows with increasing rotation (along ar- 
row 2), now at practically constant fc of ^ 350 Hz. Even 
for these models, fc is always lower than the typical peak 
frequency /max ^ 700 — 800 Hz of their spectral GW en- 
ergy densities. This is due to the specific characteristics 
of the LIGO detector, whose highest sensitivity is around 
100 Hz, thus leading to a systematic decrease of fc with 
respect to /max- 

In more rapidly rotating models, centrifugal effects 
become more important, leading to greater rotational 
deformation of the inner core, but also slowing down 
the dynamics around core bounce, ultimately limiting he 
and reducing fc (along arrow 3). Models that rotate so 
rapidly that they undergo centrifugal bounce (type IIIc) 
cluster in a separate region in the he — fe plane (along 
arrow 4), somewhat below the maximum value of he and 
at considerably lower fc- The systematics for the lower-T 
models and for other detectors is very similar. Not sur- 
prisingly, given the analogies in the two signals, a similar 
behavior of he and fc was observed in the context of ro- 
tating iron core collapse [iot - 

Figure [221 provides the same type of information shown 
in Fig. [21] but also for the advanced LIGO detector 
when the source is at 0.8 kpc {e.g., within the An- 
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FIG. 22: Location of the GW signals from core bounce in 
the hc-fc plane relative to the sensitivity curves of various 
interferometer detectors (as color-coded) for an extended set 
of models AD. The sources are at a distance of 10 kpc for 
LIGO, 0.8 Mpc for Advanced LIGO, and 5 Mpc for the Ein- 
stein Telescope. 



drom eda g alaxy) , or for the proposed Einstein Telescope 
(ET) |l26l | and a source distance of 5 Mpc. Initial LIGO 
is sensitive only to GWs coming from a moderately- 
rapidly or rapidly rotating AIC event in the Milky Way, 
but its advanced version will probably be able to reveal 
sources also outside the Galaxy, although only within the 
local group. Finally, third-generation detectors such as 
ET, may be sensitive enough to detect some AIC events 
out to ~ 5 Mpc. 



As pointed out in Sees. [Tvn and HVBl the GW sig- 
nal amplitudes and the spectral GW energy distribution 
is determined primarily by /3ic,b- Hence, given the sys- 
tematics shown in Fig. I21i one may be optimistic about 
being able to infer Pie.h to good precision from the obser- 
vation of GWs from a rotating AIC event. For example, 
as demonstrated in Fig. [THl even the knowledge of only 
/max can put some constraints on /?ic,b- However, infer- 
ring accurately the properties of the progenitor WD us- 
ing exclusively information provided by GWs may be ex- 
tremely difficult given the highly degenerate dependence 
of Ac,b on the various precollapse WD model parame- 
ters discussed in Sec. IIII Bl To elaborate on this point, 
we show in Fig. [23] the relation between /?ic,b and the 
precollapse WD parameter Pi. Even if GWs can provide 
good constraints on /3ic,b, a rather large variety of models 
with different initial rotational properties would be able 
to lead to that same /9ic,b and additional astrophysical 
information on the progenitor will be needed to deter- 
mine the precollapse rotational configuration. The only 
exception to this is the possibility of ruling out uniform 
WD progenitor rotation if /3ic,b > 18% (c/ Sec. [ml]). 
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TABLE IV: GW signal characteristics for the high-T AIC 
models: Eqw is the total GW energy, /max is the peak fre- 
quency of the GW energy spectrum, A/50 is the frequency 
interval around /max that emits 50 % of -Egw- The nonrotat- 
ing models are omitted here. 



AIC 




/max 


A/50 


model 




fHzl 


fHzl 


AUl 


1.1 


742.7 


31 


AU2 


5.7 


782.7 


28 


AU3 


7.8 


786.7 


27 


AU4 


15.7 


816.0 


49 


AU5 


17.0 


831.9 


120 


DUl 


0.1 


768.0 


343 


DU2 


0.3 


770.0 


556 


DU3 


0.4 


747.0 


473 


DU4 


0.9 


745.0 


304 


DU5 


2.0 


765.7 


26 


DU6 


4.6 
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FIG. 23: The inner core parameter Ac,b at bounce is plotted 
as a function of the precollapse parameter /3i for high-T mod- 
els. Due to different central densities and rotation profiles 
of the precollapse WD models, there is no one-to-one corre- 
spondence between /3ic,b and Pi. Hence, although one can 
extract Ac,b accurately from the bounce AIC GW signal, it 
is impossible to put strong constraints on /3i using the GW 
signal. 



PROSPECTS FOR NONAXISYMMETRIC 
ROTATIONAL INSTABILITIES 



Nonaxisymmetric rotational instabilities in PNSs 
formed in AIC or iron core collapse have long been 
proposed as strong and possibly long-lasting sources of 
GWs (see, e.g., [4l|] for a recent review). The post- 
bounce GW emission by nonaxisymmetric deformations 
of rapidly rotating PNSs could be of similar amplitude 
as the signal from core bounce and, due to its poten- 
tially much longer duration, could exceed it in emitted 
energy (e.g., 62 , [gI, [131 ) ■ Moreover, since the character- 
istic GW amplitude he scales with the square root of the 
number of cycles, the persistence of the nonaxisymmet- 
ric dynamics for many rotation periods can drastically 
increase the chances for detection. 

The simulations presented in this paper impose ax- 
isymmetry, hence we are unable to track the formation 
and evolution of rotationally induced nonaxisymmetric 
structures. Nonetheless, since the dynamical hi gh-/ ? in- 
stability can develop only at (3 above /3dyn — 0.25 [54[55t. 
we can still assess the prospects for such instabilities by 
studying the values of /3 reached by our AIC models. 
Moreover, as we shall see below, the analysis of the rota- 
tional configuration of the newly formed PNS can give a 
rough idea about the outlook also for low-/? instabilities. 

As shown in Sec. IIIIBl for not very rapidly rotat- 
ing models, the parameter /3ic,b of the inner core at 
bounce increases with the progenitor rotation and sat- 
urates at ^ 24.5% (see Fig. [TU)) . Immediately after 
bounce, the inner core re-expands and, after undergo- 
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ing several damped oscillations, settles into a new quasi- 
equilibrium state with a /3ic,pb typically smaller by ^ 3 % 
(in relative value) than that at bounce. The highest value 
of Ac,pb of our entire model set is ~ 24 % (observed in 
model AD12f4) and most other rapidly rotating mod- 
els reach values of /3ic,pb that are well below this value 
(c/. Tab. [n|. Hence, we do not expect the high-/3 in- 
stability to occur immediately after bounce in most AIC 
events. 

On the other hand, the matter around the PNS ex- 
periences rapid neutrino-cooling (not modeled by our 
approach) and the PNS contracts significantly already 
in the early postbounce phase. This results in spin-up 
and in a substantial increase of /?ic.pb- Using the VUL- 
CAN/2D code, Ott [131 studied the postbounce evolution 
of the PNS rotation of the Dessart et al. AIC models [H]. 
He found that, in the case of the rapidly rotating 1.92Af0 
model, the postbounce contraction leads to a growth of 
Ac.pb by - 50% from ~ 14 % to - 22 % in the ini- 
tial ^ 50 ms after bounce. We expect that a similar in- 
crease of the parameter /3ic,pb should take place also for 
the rapidly rotating AIC models considered here. More 
specifically, if we assume that (3 increases by ~ 50 % 
within ~ 50 ms after bounce, we surmise that AIC mod- 
els with /3ic,b ^ 17 % at bounce should reach /?ic,pb ^ Pdyn 
within this postbounce interval and thus become subject 
to the high-/? dynamical instability. 

As mentioned in Sec. IIII Bl uniformly rotating WDs 
cannot reach /3ic,pb in excess of ~ 10.5%. Hence, they 
are unlikely to become subject to the high-/? dynamical 
nonaxisymmetric instability, but may contract and spin 
up to /3 > Pscc — 14 % at which they, in principle, could 
experience a secular nonaxisymmetric instability in the 
late postbounce phase. However, other processes , e.g., 
MHD dynamos and instabilities (see, e.g., |l27l . ll28j ) may 
limit and/or decrease the PNS spin on the long timescale 
needed by a secular instability to grow. 

In addition to the prospects for the high-/? instability, 
the situation appears favorable for the low-/3 instability 
as well. The latter can occur at much lower values of 
/? as long as the PNS has si gnifi cant differential rota- 
tion (see, e.g., [stI. [ssI. [6iI. l63l. [63. Issl. [66t and references 
therein). While this instability's true nature is not yet 
understood, a necessary condition for its development 
seems to be the existence of a corotation point inside the 
star, i.e., a point where the mode pattern speed coin- 
cides with the local angular velocity [6^, Il29j . Bearing in 
mind that the lowest order unstable modes have pattern 
speeds of the order of the characteristic Keplerian angu- 
lar velocity C(richar) [111, we can easily verify whether 
such a criterion is ever satisfied in our models. Assum- 
ing a characteristic mass of the early postbounce PNS of 
~ O.8M0 and a radius of ^ 20 km, we obtain a character- 
istic Keplerian angular velocity of f2char ^ 4 rad ms~^. 
Because most AIC models that reach /?ic,pb > 15 % have 
a peak value of 17 > 5 rad ms" , it is straightforward to 
conclude that these models will have a corotation point 
and, hence, that the low-/? instability may be a generic 



feature of rapidly rotating AIC. We note that even uni- 
formly rotating precollapse models have strong differen- 
tial rotation in the postshock region outside the inner 
core. However, further investigation is needed to infer 
whether such models may be also be subject to low-/3 
dynamical instability. 

As a concluding remark we stress that the above dis- 
cussion is based on simple order-of-magnitude estimates 
and is therefore rather inaccurate. Reliable estimates 
can be made only by performing numerical simulations 
in 3D that adequately treat the postbounce delcptoniza- 
tion and contraction of the PNS and that investigate 
the dependence of the instability on /?ic,pb, on the de- 
gree of differential rotation, and on the thermodynamic 
and MHD properties of the PNS. Finally, these calcu- 
lations will also establish what is the effective long-term 
dynamics of the bar-mode deformation. In simulations of 
isolated poly tropes [si] . [ssj and from perturbative calcu- 
lations |130| , it was found that coupling among different 
modes tends to counteract the bar-mode instability on a 
dynamical timescale after its development. It is yet un- 
clear whether this behavior will be preserved also in the 
AIC scenario, where infalling material with high specific 
angular momentum may lead to significant changes. This 
will be the subject of future investigations. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have presented the first general- 
relativistic simulations of the axisymmetric AIC of mas- 
sive white dwarfs to protoneutron stars. Using the 
general-relativistic hydrodynamics code CoCoNuT, we 
performed 114 baseline model calculations, each start- 
ing from a 2D equilibrium configuration, using a finite- 
temperature microphysical EOS, and a simple, yet effec- 
tive parametrization scheme of the electron fraction Y,, 
that provides an approximate description of deleptoniza- 
tion valid in the collapse, bounce, and very early post- 
bounce phases. The precollapse structure and rotational 
configuration of WDs that experience AIC is essentially 
unconstrained. This prompted us to carry out this work. 
With our large set of model calculations, we have inves- 
tigated the effects on the AIC evolution of variations in 
precollapse central density, temperature, central angular 
velocity, differential rotation, and deleptonization in col- 
lapse. The inclusion of general relativity enabled us to 
correctly describe the AIC dynamics and our extended 
model set allowed us for the first time to study system- 
atically GW emission in the AIC context. 

We find that the overall dynamics in the collapse phase 
of AIC events is similar to what has long been estab- 
lished for rotating iron core collapse. A universal division 
in homologously collapsing inner core and supersonically 
infalling outer core obtains and the self-similarity of the 
collapse nearly completely washes out any precollapse dif- 
ferences in stellar structure in the limit of slow rotation. 
Due to the high degeneracy of the electrons in the cores of 
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AIC progenitor WDs, electron capture is pred icted to be 
strong already in early phases of collapse [2^ , leading to 
a low trapped lepton fraction and consequently small in- 
ner core masses Mic,b at bounce of around 0.3 Mq which 
decrease somewhat with increasing precollapse WD tem- 
perature due to the temperature dependent abundance 
of free protons. Test calculations motivated by potential 
systematic biases of the AIC Y^ip) trajectories obtained 
from [2^ (see Sees. Ill and Hi C 3[) with inner-core val- 
ues of Ye increased by ~ 10% and ~ 20% yielded values 
Mic,b larger by ~ 11% and - 25%. 

Our simulations show that rotation can have a pro- 
found influence on the AIC dynamics, but will always 
stay subdominant in the collapse of uniformly rotating 
WDs whose initial angular velocity is constrained by the 
Keplerian limit of surface rotation. In rapidly differen- 
tially rotating WDs, on the other hand, centrifugal sup- 
port can dominated the plunge phase of AIC and lead 
to core bounce at subnuclear densities. We find that 
the parameter /3ic,b = (^'rot/|Vr|)ic.b of the inner core at 
bounce provides a unique mapping between inner core ro- 
tation and late-time collapse and bounce dynamics, but 
the mapping between precollapse configurations and /3ic,b 
is highly degenerate, i.e., multiple, in many cases very dif- 
ferent precollapse configurations of varying initial com- 
pactness and total angular momentum, can yield prac- 
tically identical /3ic,b and corresponding collapse/bounce 
dynamics. 

Recent phenomenological work presented in [10, 113 on 
the potential EM display of an AIC event has argued for 
both uniform WD rotation [H, Il0l| and massive quasi- 
Keplerian accretion disks left behind at low latitudes af- 
ter AIC shock passage. The analysis of our extensive 
model set, on the other hand, shows that uniformly ro- 
tating WDs produce no disks at all or, in extreme cases 
that are near mass shedding at the precollapse stage, 
only very small disks (Mdisk ^ O.OSMq). Only rapidly 
differentially rotating WDs yield the large disk masses 
needed to produce the enhanced EM signature proposed 
in Plii. 

An important focus of this work has been on the GW 
signature of AIC. GWs, due to their inherently multi- 
D nature, are ideal messengers for the rotational dy- 
namics of AIC. We find that all AIC models following 
our standard Ye{p) parametrizations yield GW signals 
of a generic morph olog y which has been classified pre- 
viously as type III j60l ItsI. \7T\. This signal type is due 
primarily to the small inner core masses at bounce ob- 
tained in these models. We distinguish between three 
subtypes of AIC GW signals. Type Illa occurs for 
Ac.b ^ 0.7% (slow rotation), is due in part to early post- 
bounce prompt convection and results in peak GW am- 
plitudes |/imax| 5 X 10~^^ (at lOkpc) and emitted en- 
ergies £^GW ^ few X 10~^ Mq (? . Most of our AIC mod- 
els produce type Illb GW signals that occur for 0.7 < 
Ac.b ^ 18% (moderate/moderately rapid rotation) and 
yield 6 X 10-22 < \h^^y\ (at lOkpc) < 8x IQ-^i and emit- 
ted energies of 9 x 10-i"A/o < Eq^ < 2 x IQ-^ Mq . 



Rotation remains subdominant in type Ilia and type Illb 
models and we find that there is a monotonic and near- 
linear relationship between maximum GW amplitude and 
the rotation of the inner core which is best described by 
the power law |/imax| « lO^^^/??'^''. Furthermore, we find 
that the frequencies /max at which the GW spectral en- 
ergy densities of type Ilia and Illb models peak are in 
a rather narrow range from ^ 720 Hz to ~ 840 Hz and 
exhibit a monotonic growth from the lower to the upper 
end of this range with increasing rotation. This finding 
suggests that the GW emission in these models is driven 
by the fundamental quadrupole (^Z) mode of the inner 
core. 

In the dynamics of AIC models that reach /3ic,b ^ 18%, 
centrifugal effects become dominant and lead to core 
bounce at subnuclear densities. Such models must be 
differentially rotating at the onset of collapse and pro- 
duce type IIIc GW signals with maximum amplitudes of 
4.0 X 10^22 < |^^_^^| (atlOkpc) < 5.5 x lO-^i, emitted 
energies of V^'^^Mq (? < Eqw < 10"*^ Mq c^, and peak 
frequencies of 62 Hz < /max ^ 800 Hz. In contrast too 
type Ilia and Illb models, in type IIIc models, |/imax|, 
-EqWi and /max decrease monotonically with increasing 

Ac,b- 

Combining the information from signal morphology, 
l^maxl, Egw and /max, wc couclude that already first- 
generation interferometer GW detectors should be able 
to infer the rotation of the inner core at bounce (as mea- 
sured by /3ic.b) from a Galactic AIC event. Due to the 
degenerate dependence of /3ic,b on initial model parame- 
ters, this can put only loose constraints on the structure 
and rotational configuration of the progenitor WD. How- 
ever, the observation of an AIC with /3ic,b ^ 18% would 
rule out uniform progenitor rotation. 

Studying the configurations of the protoneutron stars 
formed in our AIC models, we find that none of them 
are likely to experience the high-/? nonaxisymmetric bar- 
mode instability at very early postbounce times. We es- 
timate, however, that all models that reach /Jjc.b ^ 17% 
will contract and reach the instability threshold within 
~ 50 ms after bounce. Less rapidly spinning models will 
require more time or will go unstable to the low-/? in- 
stability. The latter requires strong differential rotation 
which is ubiquitous in the outer PNS and in the post- 
shock region of our AIC models. AIC progenitors, due 
to their evolution through accretion or formation through 
merger, are predestined to be rapidly rotating and form 
PNSs that are likely to become subject to nonaxisym- 
metric instabilities. This is in contrast to the precollapse 
iron cores of ordinary massive stars that are expected 
to be mostly slowly-spinning objects [13, [111- We con- 
clude that the appearance of nonaxisymmetric dynamics 
driven by either the low-/? or high-/? instability and the 
resulting great enhancement of the GW signature may 
be a generic aspect of AIC and must be investigated in 
3D models. 

The comparison of the GW signals of our axisymmet- 
ric AIC models with the gravitational waveforms of the 
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iron core collapse models of Dimmelmeier et al. [49| re- 
veals that the overall characteristics of the signals are 
rather similar. It appears unlikely that AIC and iron 
core collapse could be distinguished on the basis of the 
axisymmetric parts of their GW signals alone, unless de- 
tailed knowledge of the signal time series as well as of 
source orientation and distance is available to break ob- 
servational degeneracies. 

The results of our AIC simulations presented in this pa- 
per and the conclusions that we have drawn on their basis 
demonstrate the complex and in many cases degenerate 
dependence of AIC outcomes and observational signa- 
tures on initial conditions. The observation of GWs from 
an AIC event can provide important information on the 
rotational dynamics of AIC. However, to lift degeneracies 
in model parameters and gain full insight, GW observa- 
tions must be complemented by observations of neutri- 
nos and electromagnetic waves. These multi-messenger 
observations require underpinning by comprehensive and 
robust computational models that have no symmetry 
constraints and include all the necessary physics to pre- 
dict neutrino, electromagnetic, and GW signatures. 

As a point of caution, we note that the generic type III 
GW signal morphology observed in our AIC models 
is due to the small inner-core values of Yg and conse- 
quently small inner core masses predicted by the Yf,{p) 
parametrization obtained from the approximate Newto- 
nian radiation- hydrodynamic simulations of [2^. Tests 
with artificially reduced deleptonization show that the 
signal shape becomes a mixture of type III found in our 
study and type I observed in rotating iron core collapse 
[4^ if the Ye in the inner core is larger by ~ 20%. In a 
follow-up study, we will employ Y,, (p) data from improved 
gcncral-relativistic radiation-hydrodynamics simulations 
[l02| to better constrain the present uncertainties of the 
AIC inner-core electron fraction. 

Although performed using general-relativistic hydro- 



dynamics, the calculations discussed here are limited to 
conformally-flat spacetimes and axisymmetry. We ig- 
nored postbounce deleptonization, neutrino cooling, and 
neutrino heating. We also neglected nuclear burning, 
employed only a single finite-temperature nuclear EOS, 
and were forced to impose ad-hoc initial temperature 
and electron fraction distributions onto our precollapse 
WD models in rotational equilibrium. Future studies 
must overcome the remaining limitations to build accu- 
rate models of AIC. Importantly, extensive future 3D 
radiation-hydrodynamic simulations are needed to ad- 
dress the range of possible, in many cases probably non- 
axisymmetric, postbounce evolutions of AIC and to make 
detailed predictions of their signatures in GWs, neutri- 
nos, and in the electromagnetic spectrum. 
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